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A  method  of  estimating  statistical  characteristics  of 
stable  ground  roughness  for  use  in  off- road  land  locomotion 
is  developed.  Computer  programs  for  obtaining  these  esti¬ 
mates  from  survey  data  have  been  prepared.  Several  esti¬ 
mated  power  spectral  densities  obtained  from  survey  data 
taken  at  military  installations  are  presented. 

Pertinent  feature0  of  these  estimated  power  spectral 
densities  are  internretable  in  terms  of  visual  features  of 
plotted  ground  profiles.  The  estimated  power  spectral  den¬ 
sities  may  be  approximated  by  similar  and  simple  form. 
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OBJECT 

Develop  a  statistical  method  of  estimating  characteristics  of 
stable  ground  roughness  for  use  in  off-road  land  locomotion.  Apply 
method  to  data  obtained  from  surveys. 


RESULTS 

Computer  programs  have  been  constructed  to  estimate  stable 
ground  roughness  characteristics.  Several  estimated  power  spectral 
densities  are  presented. 


CONCLUSIONS 

Power  spectral  density  methods  are  useful  for  characterizing 
quantitatively  stable  ground  roughness.  Features  of  the  estimated  p.s.d.  's 
are  interpretable  in  terms  of  visual  features  of  plotted  ground  profiles. 
Estimated  p.  s.  d.  's  may  be  approximated  by  similar  and  simple  form  in 
cases  considered. 
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GLOSSARY 


Stable  ground  roughness  is  the  variation  in  elevation  of  hard 
ground  whose  influence  on  vehicle  dynamics  remains  reasonably  constant 
over  distance,  and  is  free  of  obstacles. 

Statistical  estimation  is  the  technique  of  estimating  average  prop¬ 
erties  from  sample  data. 

Power  spectr  al  density  (p.s.d.  )  of  a  profile  measures  the  amount 
of  variation,  by  frequency  bands,  of  the  profile  height. 


Ground  profile  is  the  plot  of  survey  height  vs.  distance. 
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SECTION  I 
INTRODUCTION 

The  Midwest  Applied  Science  Corp.  was  asked  by  the  Land 
Locomotion  Laboratory,  Detroit  Arsenal,  to  identify  and  study  the 
several  areas  which  determine  the  speed  of  military  vehicles  under 
off-road  hard  ground  conditions.  Three  areas  have  been  identified 
as  pertinent  and  are  now  under  study;  they  are  specification  or  de¬ 
scription  of  stable  ground  roughness,  comfort  or  tolerance  limits 
of  humans  to  vehicle  oscillation,  and  dynamical  characteristics  of 
a  military  vehicle.  Although  results  are  usually  presented  separ¬ 
ately,  the  separate  studies  in  these  areas  must  be  synthesized  into 
vehicle  suspension  system  concepts  which  will  permit  increases  in 
vehicle  speeds. 

This  report  is  concerned  only  with  the  description  of  stable 
ground  roughness. 

The  problem  of  describing  the  roughness  of  fluctuations  in 
quantities  which  vary  with  time  or  distance  is  an  old  one.  It  has  been 
and  still  is  of  interest  in  diverse  areas  such  as  physics,  economics, 
mathematics,  engineering,  neurology,  etc.  Questions  of  the  type: 

What  is  the  composition  of  white  light?.  Will  the  stockmarket  fall  next 
week?.  Is  a  function  smooth?.  Will  noise  in  a  telephone  mask  the 
message?,  How  rough  is  the  road  a  car  can  travel  safely  at  60  m.p.h.  , 
etc.  ,  all  relate  to  the  roughness  of  quantities  which  change  with  time 
or  distance. 

Early  ways  of  describing  roughness  are  contained  in  the  motion 
of  bounded  variation,  magnitude  of  the  Fourier  coefficients  of  a  per¬ 
iodic  motion,  power  in  frequency  bands  of  light  as  determined  by  spec¬ 
trographs,  electrical  power  in  frequency  bands  of  fluctuating  electrical 
current,  etc.  Some  of  these  concepts  are  familiar  to  all  of  us. 
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Tliut  the  designer  of  military  vehicle  suspension  systems  is 
vitally  interested  in  ground  roughness,  or  more  particularly,  stable 
ground  roughness,  follows  from  the  substantial  number  of  test  courses 
at  the  various  military  installations  around  the  country  which  are  sup¬ 
posed  to  simulate  "rough  ground". 

A  military  vehicle  operating  off  roads  encounters  two  types  of 
hard  ground  roughness.  One  type  consists  of  obstacles  such  as  large 
rocks,  tree  stumps,  stream  beds,  abrupt  hills, etc.  The  other  type  con¬ 
sists  of  variations  in  elevation  which  are  stable  over  reasonably  large 
areas  and  change  but  gradually  with  distance.  The  latter  type  is  what 
we  have  called  stable  ground  roughness. 

In  general,  obstacles  require  special  driver  attention  dictated  by 
the  capabilities  of  the  vehicle.  What  is  an  obstacle  to  one  vehicle  is 
not  an  obstacle  to  another. 

Rapid  traversibility  of  the  stable  roughness  depends  upon  the  ve¬ 
hicle,  suspension  system,  and  cargo.  Again,  what  is  rough  to  one 
vehicle  may  not  be  rough  to  another  vehicle  or  to  the  same  vehicle  at 
another  speed.  Thus,  stable  ground  roughness  has  meaning  in  two  con¬ 
texts:  in  one,  only  geometry  is  involved;  in  the  other,  geometry  eoupled 
with  vehicle-driver- cargo  and  speed  is  of  importance.  Hence,  for  the 
vehicle-driver-cargo  system,  it  is  impossible  to  identify  the  roughest 
terrain  without  knowledge  of  speed  and  system  characteristics. 

Environmental  roughness  exists  and  is  now  being  studied  for  other 
than  land  vehicles.  Descriptions  of  the  roughness  of  the  sea,  atmospheric 
turbulence,  and  airport  runways  are  needed  to  predict  the  vibrational 
state  in  ships  and  airplanes.  Starts  have  been  made  in  obtaining  descrip¬ 
tions  in  terms  of  power  spectral  density  and  some  of  the  preliminary  re¬ 
sults  obtained  from  such  descriptions  appear  useful  to  the  designer.  We 
note,  of  course,  the  lack  of  obstacles  in  air  and  sea  travel,  and  the  lack 
of  long  unimportant  trends  such  as  gradual  hills.  Here,  as  with  land 
vehicles,  what  is  rough  to  an  airplane  or  ship  depends  upon  speed  and 
dynamical  characteristics. 
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It  is  worth  explicitly  noting  that  there  are  no  prior  references 
to  be  cited  in  connection  with  statistical  descriptions  of  stable  off¬ 
road  hard  ground  roughness. 
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OBJECT 

The  object  of  this  report  is  to  present  statistical  methods  of 
describing  or  estimating  quantitatively  the  characteristics  of  stable 
roughness  of  hard  ground  of  interest  in  vehicle  dynamics,  and  to 
present  estimates  of  roughness  obtained  from  survey  data.  Smooth¬ 
ing  techniques  must  be  developed  to  eliminate  features  of  roughness 
unimportant  in  vehicle  dynamics,  such  as  long  gradual  hills.  Further, 
the  assumption  that  the  underlying  roughness  is  reasonably  stable 
over  lengths  large  in  comparison  with  vehicle  lengths  must  be  verified. 
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RESULTS 

Figures  1-4  present  line  profiles  obtained  from  elevation  sur¬ 
veys  taken  at  the  indicated  places.  Horizontal  distance  was  measured 
with  a  tape;  elevation  was  measured  with  level  and  rod. 

In  each  of  the  first  three  figures,  the  complete  profile  is  shown; 
each  data  point  (represented  by  a  dot)  is  shown,  and  the  data  points 
are  two  feet  apart.  Figure  4  shows  every  tenth  data  point  and  the  data 
points  are  one  foot  apart.  The  vertical  scale  is  substantially  different 
from  the  horizontal;  a  twenty  foot  interval  is  marked  to  give  an  idea  of 
the  horizontal  scale.  In  Figure  4,  a  complete  section  has  been  added  at 
each  end  to  indicate  details. 

The  effect  of  ground  roughness  upon  a  vehicle  is  dependent  on  the 
vehicle's  size  and  speed.  It  is  obvious  that  in  the  context  of  today's 
military  vehicles  some  of  the  variations  in  ground  elevation  in  these 
figures  are  of  such  a  scale  as  not  to  be  of  interest  in  vehiele  dynamics. 
For  example,  in  Figure  4,  the  profile  has  an  upward  sloping  trends 
whieh  indicates  the  aseent  of  a  hill.  Superimposed  upon  this  hill  are 
the  variations  in  elevation  or  roughness  which  a  vehicle  will  feel.  In 
the  other  figures,  such  a  striking  trend  is  not  present  although  there 
are  large  scale  elevation  variations  whieh  might  be  termed  hills  (in 
these  figures)  and  are  not  of  interest  in  this  study.  However,  if  one 
mentally  eliminates  these  slowly  varying  elevation  changes,  one  notes 
that  the  remaining  variations  or  roughness  patterns  appear  to  be  fairly 
constant.  What  is  needed  is  a  quantitative  description  of  the  roughness 
pattern. 

Whatever  the  description  adopted  for  roughness  measurement, 
it  is  essential  that  it  be  usable  by  the  suspension  system  designer.  Moxe- 
over,  it  must  classify  as  having  the  same  or  different  roughness  (two 
roads  or  fields  classified  by  experience  as  having  the  same  or  different 
roughnesses).  We  shall,  in  this  report,  be  concerned  with  the  power 
speetral  density  description  of  ground  roughness. 
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A  complete  description  of  power  spectral  density  analysis 
of  roughness  on  a  line  is  given  below  for  the  interested  reader. 

Briefly  and  loosely  speaking,  the  power  spectral  density  of  an  os¬ 
cillating  function  or  ground  elevation  along  a  line  provides  a  measure 
of  the  power  or  intensity  in  frequency  (cycles /ft.  )  bands;  all  os¬ 
cillating  functions  having  the  same  p.  s.  d.  are  classified  as  having 
the  same  roughness. 

Figures  5  to  8  present  the  p.  s.  d.  ’s  (power  spectral  densities) 

estimated  from  the  profiles  given  in  Figures  1  to  4,  respectively. 

3 

The  ordinate  has  units  of  feet  and  the  abscissa  is  frequency  with  unit 
feet  1. 

The  general  features  of  these  p.  s.d.  's  are  mueh  the  same. 
There  is  a  sharp  spike  at  the  origin  with  a  more  or  less  rapid  drop 
as  frequency  increases,  ana  a  long  low  tail.  In  two  eases,  Figure  5 
and  Figure  6,  there  is  a  bump  in  the  rapid  drop-off  at  approximately 
X  =  .  025  and  X  -  .  035,  respectively.  The  tails  show  irregularities  at 
a  very  low  level,  in  all  cases  below  0.  30.  Some  of  the  irregularities 
m  the  tails  are  linked  to  inconsistencies  in  the  profile  data.  However, 
the  main  features  of  interest  in  vehiele  dynamies  lie  in  the  frequency 
range  where  the  p.  s.d.  is  descending  to  the  tails. 

Let  us  return  to  Figure  1.  (Aberdeen)  This  profile  shows  no 
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prominent  hills.  All  elevations  lie  between  -  4  feet.  This  is  more  or 
less  characteristic  of  the  area  within  the  Perryman  Mud  Course  where 
the  profile  was  taken.  Qualitatively,  the  roughness  seems  to  eonsist 
of  changes  in  elevation  of  approximately  2  feet  within  distances  of  20 
to  30  feet.  In  the  first  third,  there  is  a  fairly  periodic  oseillation  of 
ground  elevation  with  a  period  of  30  to  40  feet;  in  the  middle  third, 
they  arc  almost  gone.  One  notices  occasionally  in  the  profile  sharp 
irregularities  in  the  elevation  between  adjacent  neighboring  points; 
these  are  due  to  ruts,  holes,  rocks,  errors  in  data  recording,  etc. 
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The  Aberdeen  p.  s.d.  ,  Figure  5,  indicates  just  a  hint  of  a  bump  at 
X  -  .  025,  corresponding  to  the  more  or  less  regular  oscillation  of 
approximately  40  ft.  wave  length  in  the  profile  noted  earlier.  We 
see  below  that  the  p.  s.d.  of  the  first  half  has  a  much  more  pro¬ 
nounced  bump  at  X  =  0.025  while  the  second  half  does  not  have  this 
bump. 

The  profile  from  Ft.  Knox,  Figure  2,  also  does  not  exhibit 
any  substantial  long  range  variation,  except  after  1600  feet.  As  the 
profile  was  taken  in  bottom  land  bordering  the  Ohio  River,  this  is 
reasonable.  The  roughness  has  amplitudes  which  are  about  half  that 
observed  in  the  Aberdeen  profile.  There  is  a  pronounced  oscillation 
between  400  ft.  and  1600  ft.  with  a  wave  length  of  approximately  27 
feet;  this  explains  the  presence  of  the  bump  in  the  corresponding 
p.s.d.  estimate.  Figure  6,  at  X  =  0.035.  The  minor  irregularities 
in  file  profile  are  distributed  through  it  as  in  the  Aberdeen  profile. 

The  profile  taken  at  Yuma,  Figure  3,  exhibits  major  long  range 
changes  in  elevation,  varying  between  +12  feet  and  -10  feet.  These 
are  represented  by  hills  of  approximately  400  feet  length;  throughout 
the  profile  there  are  smaller  hills  of  about  50  feet  length;  the  perti¬ 
nent  roughness,  from  the  point  of  view  of  vehicle  dynamics,  exhibits 
no  significiant  periodicities.  The  scale  of  this  latter  roughness  appears 
smaller  than  present  in  either  Aberdeen  or  Knox.  The  p.  s.d.  estimate, 
in  Figure  7,  indicates  a  much  higher  value  near  the  origin,  correspond¬ 
ing  to  tlie  smaller  hills,  and  a  sharper  drop  due  to  less  roughness. 

The  ordinate  of  the  p.s.d.  atX  =  0.03  is  1.0  whereas  Aberdeen  and 
Knox  are  substantially  higher,  confirming  this  observation. 

All  6201  data  points  of  the  Battlefield  Day  profile  (taken  at  Fort 
Knox)  could  not  be  shown  on  the  main  profile,  the  spacing  between 
points  shown  is  10  ft.  ;  consequently  the  relevant  roughness  is  largely 
masked.  Tv\o  sections  at  either  end  have  been  plotted  completely  to 
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give  samples  tf  the  type  of  roughness  met.  The  main  profile  exhibits 
one  long  rising  trend  with  superimposed  minor  hills.  The  scale  of  the 
roughness  of  interest,  indicated  only  in  part  in  the  detailed  portions, 
appears  to  be  smaller  even  than  in  Yuma.  Rarely  do  these  roughness 
variations  exceed  1/4  feet.  The  p.s.d.  ,  Figure  8,  drops  faster  than 
in  the  Yuma  p.  s.d.,  as  would  be  expected.  The  fluctuation  in  the  tail 
of  this  p.s,.  d. ,  as  well  as  its  elevation,  will  be  discussed  below;  essen¬ 
tially  it  is  traceable  to  inconsistencies  in  the  data. 
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CONCLUSIONS 

Our  main  conclusion  is  that  the  principal  problems  of  charac¬ 
terizing  ground  roughness  by  p.s.d.  methods  are  behind  us.  We  must 
now  turn  to  the  more  difficult  task  of  relating  the  measurements  to 
the  suspension  system  design. 

The  raw  data  for  the  one  and  two  dimensional  cases  are  available 
in  IBM  card  decks.  Anyone  interested  in  these  data  may  obtain  them  by 
writing  to  MASC.  These  data  will  also  be  tabulated  and  presented  in  a 
later  report. 

The  following  practical  conclusions  appear  to  be  consistent  with 
the  results  obtained  to  date: 

One -dimensional 

a.  Visual  observed  roughness  characteristics  of  the 
profile  are  related  to  characteristics  of  the  esti¬ 
mated  p. s.d. 

b.  The  estimated  p.  s.d.  's  may  be  approximated  with 
reasonable  accuracy  by  a  relatively  simple  class 
of  functions. 

c.  Abnormal  height  and  periodic  oscillations  in  the 
tails  of  the  p.s.d.  estimates  are  traceable  to  in¬ 
consistencies  in  the  survey  data. 

d.  Sections  of  ground  which  appear  visually  to  have 
constant  roughness  characteristics  satisfy  reason¬ 
ably  well  the  stationarity  assumption  on  intervals 
of  1000-2000  feet  length. 

e.  The  normality  assumption  for  the  smoothed  data  is 
reasonably  well  satisfied  although  there  is  a  con¬ 
sistent  deviation  from  it. 
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f.  There  may  be  considerable  difficulty  in  removing 
unwanted  hills  (long  range  trends)  from  the  data. 

Running  average  methods  have  proved  satisfactory 
with  present  data. 

g.  The  one-dimensional  p.s.d.  estimate  computer  program 
is  almost  ready  for  distribution. 

Two-dimensional 

h.  P.  s.d.  estimates  obtained  to  date,  although  not  com¬ 
pletely  satisfactory,  display  general  features  ob¬ 
served  in  one- dimensional  results. 

In  general,  successes  to  date  in  obtaining  p.  s.d.  estimates  are 
quite  gratifying.  There  are  details,  however,  mostly  related  to  smoothing 
and  interpretation  for  two-dimensional  p.  s.d.'s,  which  must  be  worked 
out.  We  anticipate  no  substantial  difficulties  in  cleaning  up  these  details 
and  arriving  at  programs  available  for  distribution. 
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RECOMMENDATIONS 

•  More  surveys  are  necessary  to  provide  p.s.d.  estimates 
for  a  wider  variety  of  ground  types. 

Smoothing  procedures  must  be  examined  for  their  effects  on 
characteristics  relevant  to  vehicle  dynamics. 

Cross-correlation  effects  from  the  two-dimensional  data 


must  be  ascertained. 
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Section  II 

One -dimensional  Power  Spectral  Densities 

1.  Introduction,  Interpretation  of  Spectra 

The  spectral  analysis  of  random  functions  of  one  variable 
is  the  subject  of  a  number  of  books  and  articles  in  the  engineering 
literature.  We  will  review  in  this  chapter  the  principal  ideas  with 
emphasis  on  the  descriptive  aspects. 

It  is  often  the  case  that  a  random  function  is  of  interest  as 
an  input  to  a  linear  or  almost  linear  system.  Its  effect  in  such  a 
problem  is  conveniently  separated  into  the  sum  of  effects  due  to  fre¬ 
quency  components.  The  effect  of  terrain  roughness  on  the  vibration 
of  a  traveling  vehicle  is  such  a  problem.  As  a  first  approach  we  will 
consider  a  single  track  vehicle.  The  profile  heights  on  the  path  of 
the  vehicle  will  be  a  function,  h(x),  of  the  distance,  x  along  the 
track.  The  roughness  of  the  profile  will  be  a  matter  of  the  variation 
of  h(x)  about  some  smooth  medial  terrain  consisting  of  the  large 
hills  etc.  The  problem  of  obtaining  this  smooth  terrain  is  a  practi¬ 
cal  one  and  is  discussed  in  a  later  section.  Let  us  denote  the  devia¬ 
tions  of  the  profile  height  by  d(x).  These  will  average  zero  and  their 
variation  is  what  we  mean  by  roughness.  As  a  drastic  simplification 
of  separation  into  frequency  components,  we  might  have 

M  m 

d(x)=  £  (  ak  sina>kx  +  bRCOS  u>kx  )  • 

K=0 

(l.D 

It  is  very  unrealistic  to  consider,  as  in  (1.  1),  that  only  a  few  frequen¬ 
cies  are  present,  we  should  admit  all  frequencies  each  in  an  infinitesi¬ 
mal  amount.  These  ideas  are  familiar  to  most  engineers  as  an  intro- 


*Generally  we  will  distinguish  between  frequencies  in  radians  per  foot 
(u)  and  cycles  per  foot  (X). 
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duction  to  the  Fourier  integral.  Here,  however,  the  fact  of  ran¬ 
domness  complicates  matters.  If  d(x)  with  fixed  frequencies  is 
random,  then  it  is  the  coefficients  a^.b^  that  are  random.  How  to 
deal  with  infinitesimal  random  quantities  in  the  case  of  continuous 
frequencies,  presents  a  special  problem  in  probability  theory  which 
we  cannot  hope  to  deal  with  here.  The  unrealistic,  fixed  frequency, 
case,  does  "converge  to"  the  realistic  case  and  it  is  worthwhile  to 
consider  it  in  some  detail. 

If  a  function  of  the  form  (1.  1)  is  used  as  an  input  to  a  linear 
system,  the  output  also  has  this  form 


a(x)  =  Z  { sin  cuk  x  +  b£  cos  x  ) 


(1.2) 


The  coefficients  a^.b^  of  the  k'th  pair  depend  only  on  a^b^.  of  the 
k'th  pair  of  (1.  1).  The  form  of  dependence  is  stated  most  clearly 
if  we  change  the  form  of  (1. 1),  (1.2) 

an  sin  (dn  x  +  bn  cos  ojn  x  =  An  cos  (a»n  x  +  <£n) 

(1.3) 

An  8  +  bn  c£n  =  arctan  (an/bn) 


Then  with  the  general  term  of  (1.2)  in  the  same  form  it  may  shown 
that 


An  =  An  P  (  CUn) 

4>'n  =  <£n  +  $  (tU"> 


(1.4) 


The  functions  P  (to)  and  $  (u)  combine  to  form  what  is  called  the 
transfer  function  of  the  linear  system.  The  phase  of  the  input  does 
not  affect  the  amplitude  of  the  output,  nor  does  the  amplitude  of  the 
input  affect  the  phase  of  the  output.  In  many  problems  the  ultimate 
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effect  does  not  depend  on  the  phase  of  the  output  so  that  only  A'  and 
hence  only  An  is  of  interest.  The  various  An  arranged  according 
to  the  sizes  of  the  u>n  form  what  is  called  the  frequency  spectrum  of 
the  input.  Their  sizes  give  the  amounts  of  the  frequencies  present. 

For  a  continuum  of  frequencies  this  concept  becomes  that  of  a  spectral 
density. 

This  is  a  quick  account  of  the  standard  method  of  analyzing 
the  effect  of  an  input  d(x)  on  a  linear  system.  We  have  not  mentioned 
that  d(x)  is  random  --by  random  we  mean  that  d(x)  may  turn  out  to 
be  any  one  of  an  "ensemble"  of  functions  and  we  have  no  hope  of  deter¬ 
mining  which  one,  but  we  can  determine  or  estimate  some  ensemble 
averages.  We  will  assume  that  only  the  coefficients  an>  are  random 
in  the  expression  (1.1).  This  is  no  restriction,  for  we  can  assume 
that  all  the  possible  frequencies  are  represented  among  the  u' s  and 
that  their  absence  or  presence  is  a  matter  of  which  coefficients  are 
or  are  not  zero. 

It  will  frequently  be  convenient  to  use  the  exponential  form  of 
the  trigonometric  series  (1.  1).  If  the  uk  are  given  for  k=0,  1 , 2,  .  .  .  ,  M 
and  uq  =  )  this  is 

d(x)--£  ak^x 
K=-M 

where 

CZ|,=  4'(b|<— iok)  k>0 

(1.5) 

a_k  =  CL*  =  y  (  bk  +  iak)  k  >  0 

cl  o  -  b0 

Complex  quantities  appear  on  the  right,  but  d(x)  is  real. 
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We  assume  two  main  properties  for  the  randomness  of 
d(x).  The  first  is  that  the  ensemble  average  at  any  point  x  is 
zero.  From  the  linear  property  of  the  expected  value  (ensemble 
average)  operation, 

E  (d  (x)}=  £  E  {  aj  e'Wk*  =  0  (1.6) 

The  series  above,  being  identically  zero  as  a  function  of  x,  must 
have  zero  coefficients,  i.e. 

E  {“4  =T<E{»k}±  iE{ok}=0 

E  {°k}  •  E  M  ■  0 


Our  second  assumption  is  that  of  stationarity.  A  consequence  of 
the  assumption  is  that  any  expected  value  concerning  the  d's  must 
remain  unchanged  for  different  placements  of  the  origin  of  x.  In 
particular  the  following  does  not  depend  on  x: 


E  {d(x  +  s)  d(x)  }  =  E  {  X  laka*eiw^*5,'iw)’ 


-  I  l  E 
k  J 


i(cok-u>i)x  iiokS 


{ak  a,* 


In  this  expression  we  have  made  use  of  the  fact  that  d(x)  =  d*(x) 
for  convenience.  The  series  on  the  right,  which  must  be  identi¬ 
cally  constant  for  all  x,  falls  into  two  parts, 


x 

k 


E  {  ak  a*  }  eiw^X.E{Mi*}ei(Vi,se 

k^J 


,u>ks 
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If  the  Ujj  are  all  different,  which  we  may  assume  from  the  start, 
then  the  second  series  must  be  identically  constant  as  a  function 
of  x  hence  its  coefficients  must  be  zero; 

E  {ak  Oj*  }  -  0  k  ) 

Thus  if  d(x)  is  to  be  stationary,  it  is  necessary  that  the  complex 
coefficients  be  mutually  uncorrelated.  In  this  case 

E  {d(x  +  s)  d(x)}  =  R(s)=£E{aka*}eiwk8  (1.7) 

K 

R(s)  is  called  the  covariance  function  of  d(x).  From  (1.3) 

2 

laJ  '  ak  ak  ■  T  <°J  +  bk  >  ■  i  aJ  (taking  A.k*  Ak) 

laj2*  ■  A5 


inis  gives 


y;i«V  M  f  , 


(1.8) 


which  links  the  covariance  fuction  with  the  spectrum.  Note  that 
the  components  are  all  in  phase,  in  fact  we  can  write 


R(s)=  E  {a*}  +  \  I  E  {a^}  cos  o;ks 


(1.8’) 


Therefore  R(O)  is  a  maximum  (positive)  value  of  R(s).  The  co¬ 
efficients  in  the  trigonometric  series  (1.8)  are  ensemble  averages 
and  are  possible  to  know  or  to  estimate.  When  they  are  arranged 
in  order  of  the  sizes  of  the  u^,  they  form  what  is  called  the  power 
spectrum  of  the  random  function  d(x).  From  (1.8)  The  power 
spectrum  is  the  frequency  spectrum  of  the  covariance  function. 
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To  summarize  what  we  have  discussed  to  now,  we  will 
deal  with  random  inputs  to  linear  systems.  The  spectrum  will 
characterize  the  effect.  However,  the  spectrum  is  random  and 
cannot  be  known  except  for  its  ensemble  averages.  The  power 
spectrum  is  a  collection  of  such  averages.  The  different  spectral 
values  are  uncorrelated.  The  power  spectrum  is  the  spectrum 
of  the  covariance  function. 

To  interpret  the  power  spectrum,  let 
frequency  values  are  very  close  together, 
that 

R(0)  =  E  {A5}+  !  E  {4- 

From  (1.7)  K_1 

R(0)  =  E  {d2(  x)} 

Assuming  that  d(x)  is  approximately  normal  and  has  a  mean  of 
zero.  R(O)  is  its  variance  and  the  probability  is  about  .  003  that 
d(x)  will  exceed  in  absolute  value  3\/R(0).  Thus  R(O)  is  a  meas¬ 
ure  of  the  approximate  extremes  of  sizes  in  profile  height.  The 
usual  values  of  d(x)  lie  between  +  From  (1.9)  we  see  that 

the  left  side  is  limited  by  the  range  of  oscillation  of  d(x),  so  that  the 
terms  of  the  sum  of  the  right  side  are  also  limited  in  size.  If  there 

2 

are  many  frequencies,  then  most  of  them  have  small  values  for  E  (A  ^}. 

Equation  (1.  10)  has  another  interpretation  through  the 
ergodic  hypothesis  which  says  that  ensemble  averages  may  be  found 
as  limits  of  time  averages.  In  particular 


us  imagine  that  the 
We  have  from  (1.  8) 


Ai} 


(1.9) 


(1.  10) 


Lim  -f  J^T(d  ( x  ))Zdx  =  E  {(d(x))2}  =  R  (0) 
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The  left  side  is  a  measure  of  the  average  variability  of  d(x)  over 
a  long  stretch.  By  equation  (1. 9)  this  is  divided  into  averages 
of  uncorrelated  components  each  of  which  has  a  separate  effect  on 
the  linear  system.  In  our  work  we  will  usually  be  interested  more 
in  frequency  bands,  and  it  is  natural  to  consider  that  the  total  vari¬ 
ance  can  be  partitioned  into  partial  variances  due  to  frequency  bands. 
In  fact,  suppose  we  have  several  random  functions  d^x),  d2(x), 
etc.  each  of  which  may  be  represented  as  a  trigonometric  series  of 
form  (1. 5)  with  random  coefficients.  Since  the  frequencies  are  dif¬ 
ferent  we  will  assume  the  coefficients  to  be  uncorrelated.  Thus 
d^(x),  d„(x),  etc.  are  also  uncorrelated  so  that  the  variance  of  the 
sum  is  the  sum  of  the  variances 

Vor  {  d,  (x)  +d2(x)  +  ...  }=  R|  (o)  +  r2  (0)+  ... 

Calling  the  sum  d(x),  we  have 

Vor{dM}=£iE{4.A5}+I2E{.^}  +  ...  (1.12) 


which  shows  that  the  variance  of  d(x)  can  be  partitioned  in  to  the 
sums  of  variances  of  various  frequency  components  each  of  which 
may  be  thought  of  as  a  separate  random  function  the  sum  of  which  is 
d(x)  itself. 

Sometimes  a  power  spectrum  has  a  bump  which  we  can 
"isolate"  by  sketching  a  "normal"  spectrum  on  which  it  is  superim¬ 
posed.  The  power  in  the  bump  is  the  added  area  over  the  normal 


*  Sj,  for  example,  means  sum  over  those  u^'s  which  are  frequency 
components  of  d^(x),  etc.  . 
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spectrum.  We  may  visualize  this  power  as  being  manifested  ap¬ 
proximately  by  a  sine  wave  of  the  frequency  at  the  eenter  of  the 
bUmp  and  amplitudes  which  may  change  but  whieh  average  about 
the  square  root  of  twice  the  area.  With  pronouneed  bumps  such 
oscillations  are  often  visible  in  the  profile  data. 


2.  Relation  of  the  theory  to  the  estimation  problem. 

In  estimating  speetral  values  from  profile  data  we  will 
be  motivated  principally  by  two  relations.  The  first  is  (1.8) 
which  links  the  power  speetrum  to  the  eovarianee  function,  and 
the  seeond  is  another  aspect  of  the  ergodie  hypothesis  (1.  11), 

R(s)  =  Lim  y  f  fd(x  +  s)  d  ( x )  1  dx  (2.1) 

T  *o  L  J 

These  equations  suggest  that  the  values  of  R(s)  might  be  esti¬ 
mated  by  taking  mean  lagged  produets  of  the  data,  and  then  the 
values  of  the  power  speetrum  might  be  gotten  by  determining  the 
Fourier  coefficients  in  the  expansion  of  R(s).  Before  going  into 
details  we  must  make  some  things  definite  that  have,  to  now,  been 
left  vague. 

We  will  assume  a  continuous  power  speetrum  given  by  a 


spectral  density  function  so  that  (1.8)  beeomes 
r00  .  2  7TiXs  ,  r 00 


.  2  tT\  as  r  u-! 

f(A)e  dX  =  I  f  (X)  cos  27r\s  d  X 


( f  (-X)  =  f  (X)) 


(2.2) 


The  quantity  f(ujc)dX  may  be  associated  with  the  value  of  the  sum 
1  2 

of  An}  for  those  un  within  a  band  of  width  dX  about  Uj. 

The  dimension  of  the  values  of  f(X)  is  (distanee)^  sinee  dX  is  in  eyeles 
per  foot  or  (dist.  )  '  and  is  (dist.  )^. 
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We  must  also  take  into  account  the  form  of  the  data  that 
we  will  have  to  work  with.  This  will  be  n  values  of  an  outcome 
of  d(x)  taken  at  evenly  spaced  intervals: 

d  ( x0  +  k  A)  ,  k=0,l,2,...,n-l. 

We  shall  define  a  mean  lagged  product  as 

n-j-l 

X  d(x0  +  kA+ jA)d(x0  +  kA)  if  j>0 

k  =  0 
n-l 

£  d(x0+ kA  +  jA)d(x0+kA)  if  j<0 

k  =  j 


the  definition  of  R(s),  (1.7)  that 

E  {r,-  }  =  R  (i A) 

By  the  ergodic  hypothesis,  the  probability  limit,  as  n°o  of  r^ 
is  R(jA).  It  is  also  evident  that  we  cannot  hope  to  estimate  any 
values  of  the  covariance  function  other  than  R(jA)  and  that  this 
limitation  is  imposed  by  the  fact  that  we  have  only  observed  evenly 
spaced  data.  Our  data  is  also  limited  to  a  finite  number  of  values 
so  that  the  number  of  covariance  values  is  also  limited.  We  should 
limit  the  number  of  lags  so  that  the  number  in  terms  (n-  |j|)  in  the 
mean  (2.  3)  is  sufficient  to  make  it  a  good  estimate.  Let  us  denote 
the  maximum  lag  by  mim  is  calculated  for  j  =  0,  +  1,  +  2.  .  . +m. 

The  above  indicates  how  we  shall  estimate  the  covariance 
function.  The  determination  of  its  Fourier  coefficients  is  accom¬ 
plished  by  taking  a  linear  combination  of  its  values  with  cosine  co- 


n-j 

I 


then 


It  is  apparent  from 
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efficients.  We  will  investigate  the  properties  of  linear  combinations 
of  the  rj's  with  any  sort  of  coefficients  Wj 

m  rn 

W  *  Z  ri  wj  (2.6) 

i  =  -m 

The  expected  value  of  this  estimate  is 

E  (w  }  =  Z  R(jA)wj 

^  J  -m 

=  ^ /”f(X)wicos  2  7rXjAdX 

-m  -co  1 

Wj  cos 27rXj  A  jdX  (2.7) 

*  /  f  (X)W(X)dX 

where 

W(X)  =  Z  wi  cos27rXj  A  (2.8) 

-  m 

Ideally  if  in  (2 .  .7 )  E  {^}  were  to  be  the  value  f(  XG)  then  W(X) 
should  be  a  Dirac  impulse  or  "delta"  function  centered  at  XG.  How¬ 
ever  W(X)  is  a  finite  Fourier  series  and  cannot  be  a  Dirac  Function 
--in  fact  it  is  periodic  with  a  period  of  1/A  and  cannot  be  zero  ex¬ 
cept  at  isolated  points. 

The  function  W(X)  of  (2.8)  is  called  the  "spectral  window' 
The  estimates  $7  obtained  by  linear  combinations  of  the  covariance 
estimates  are  not  estimates  of  the  values  of  the  power  spectrum, 
but  of  averages  of  those  values  of  f(X)  admitted  by  the  window  as 
in  (2.7).  By  making  W(X)  have  a  narrow  peak  of  unit  area  near 
X0  and  very  low  values  elsewhere  we  will  estimate  very  nearly  the 
power  over  an  interval  of  frequencies  near  XQ.  This  corresponds 
to  one  of  the  sums  in  the  partition  (1. 12)  of  the  total  power. 

The  properties  of  W(X)  arising  from  its  form  as  a  finite 
Fourier  series  (2.8)  give  several  limitations  of  the  linear  estimate 
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From  the  fact  that  W(X)  is  even  and  periodic,  W(X  )  =  W(l/A  -  X  ). 

o  o 

This  means  that  while  trying  to  estimate  f (X0)  we  will  unavoidably 
pick  up  what  power  there  is  in  f(X)  around  l/'A  -  XQ  since  there  is 
also  a  peak  in  the  window  around  this  point.  All  the  power  beyond 
the  frequency  l/2A  must  unavoidably  be  confused  with  power  in  the 
range  0,  l/2A.  There  is  also  another  peak  at  -X  which  picks  up 
the  power  around  f(-Xo),  but  since  the  spectral  density  function  is 
even,  this  is  not  a  bad  feature;  it  means,  however,  the  E  {  W  }  ap¬ 
proximates  f(XQ)  +  f(-X  )  unless  the  total  area  of  the  two  peaks  is 
unity.  The  window  has  an  infinity  of  other  peaks,  all  of  them  con¬ 
tributing  spurious  power  to  the  estimate,  in  the  same  manner  as  the 
peak  at  l/A  -  Xq.  These  spurious  contributions  are  called  "aliasing", 
and  their  exact  nature  depends  on  the  sampling  interval  A.  In  de¬ 
signing  the  sampling  method,  A  must  be  chosen  so  that  the  total 
power  in  the  frequencies  beyond  1/2A  is  not  appreciable. 

From  the  fact  that  W(X)  can  be  zero  only  at  isolated  points 

m  1 

(actually  z  w(  2  vriA  log  z)  is  a  polynomial  in  z  of  degree 
2m  +  1  and  cannot  have  more  than  2m  +  1  zeros)  it  follows  that  in 
trying  to  estimate  f(XQ)  we  will  unavoidably  pick  up  some  power 
in  f(X)  near  the  non-zero  parts  of  W(X).  If  the  spectrum  has  a 
high  spot  at  one  of  the  off  center  non-zero  places  of  W(X)  the  power- 
picked  up  may  be  considerable .  If  W(X)  is  negative,  this  spurious 
power  is  subtracted  from  the  total  and  may  make  the  expected  spectral 
estimate  go  negative.  These  non-zero  off  center  parts  of  the  window 
are  called  "side  lobes".  To  make  them  small  one  must  at  the  same 
time  make  the  peak  broader  so  that  E  {W}  becomes  an  average  of  a 
wider  band  of  frequencies.  The  band  width  is  inversely  proportional 
to  m  and,  though  the  exact  definition  of  band  width  is  not  fixed  (the 
window  cannot  be  a  rectangle)  it  is  usually  taken  to  be  about  l/2mA. 

In  this  way  the  size. of  the  sample  is  related,  through  the  maximum 
lag  number,  m  ,  to  the  resolution  or  band  width  of  the  estimate. 
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These  limitations  on  the  exactness  of  the  estimate  are 
not  of  a  statistical  nature.  The  statistical  errors  are  random  de¬ 
viations  of  the  estimate  W  from  its  expected  value  E{W}  (2.6), 
(2.7).  What  we  have  considered  so  far  are  properties  of  E{W}; 
it  is  an  average  of  spectral  power  over  a  band  of  frequencies  with 
unavoidable  contamination  from  aliasing  and  from  the  side  lobes  of 
the  window.  Considerable  attention  has  been  paid  to  the  shape  of 
the  window  (see,  for  example,  Blackman  and  Tukey,  "The  Measure¬ 
ment  of  the  Power  Spectra",  Dover  1958,  Parzen,  "Mathematical 
Considerations  in  the  Estimation  of  Spectra",  Technometrics  Vol.  3, 
No.  2,  May  1961,  167-190).  Most  of  the  considerations  in  the  se¬ 
lection  of  a  window  are  more  a  matter  of  art  than  of  hard-eyed  analy¬ 
sis,  but  some  of  the  aspects  of  calculation  are  based  on  simple 
mathematical  identities  which  we  will  now  present. 

The  first  window  that  comes  to  mind  is  the  finite  Fourier 
series  whose  coefficients  are  the  first  m  Fourier  coefficients  of 
the  Dirac  function  centered  at  X  .  In  terms  of  (2.8) 


which  gives  the  window  function 
_  m 

W  (X)  =  A  Y.  cos  27™oJ  A  COS  27TAj  A  (2.10) 

-m 


This  sum  may  be  expressed  in  closed  form  by  writing  the  cosines 
in  exponential  form  and  summing  the  resulting  geometric  series, 
or  by  using  the  formula  2  cos  j©  =  s^n  (m+l/2)  0/sin  G/2: 


-  m 
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W  (X)  -  JL  l  sin  tt{  2m  +  l  )(X+X0)  A  sin77~(2m+l)  ( X-  X0)  A 
2  \  sin7r(X  +  X0)A  sin7T(X-X0)A 

(2.1 


This  function  has  high,  sharp  peaks  at  Aq  and  -A  as  well  as 
the  aliasing  peaks  at  -A  +  k/ A.  It  also  has  large  side  lobes 
including  negative  ones  (Graphs  may  be  seen  in  Blackman  and  Tukey, 
referred  to  above,  or  in  Jenkins,  "General  Considerations  in  the 
Analysis  of  Spectra,  Technometrics  Vol.  3,  No.  2,  May  1961,  133- 
166).  Both  of  these  are  bad  features  of  this  window.  --We  have  al¬ 
ready  diseussed  the  side  lobes,  the  narrowness  has  the  effect  of 
increasing  the  statistical  errors  as  we  shall  see  later.  Integrating 
(2J0)  term  by  term  over  (-  1/2A,  1/2A),  all  but  the  0  term  vanishes 
and  the  result  is  unity.  This  means  that  the  two  peaks  of  w  in 
this  interval  have  a  total  area  of  unity  and  for  A  =  0  the  single  peak 
has  unit  area.  Thus  w  always  estimates  an  approximation  to 


f(xo). 

All  of  the  spectral  windows  presented  in  the  literature  are 

modifications  of  this  simple  window- -actually  of  these  simple  win- 

dows  sinee  there  are  different  ones  for  each  A  q.  The  modification 

does  not  depend  on  A  and  consists  in  multiplying  the  coefficients 

w.  of  the  simple  window  (2.9)  by  constants  c.  ,  (c  .=c.).  These 
J  J  ~  J  J 

may  be  applied  to  the  covariances  m  before  the  calculations  (2.6) 

are  made  with  the  simple  window  coefficients  so  that  the  simple 

windows  for  each  Aq  are  all  applied  to  the  same  modified  covariances. 

The  modified  window  function  may  be  calculated  in  terms  of  the  simple 

window.  We  will  use  the  notation  w  and  Vv  to  refer  only  to  the 

simple  window  defined  by  (2.  9)  and  (2.  10).  We  will  express  (2.  10) 


in  exponential  form: 
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Let  us  denote  the  cosine  transform  of  the  constants  c^  by 

m  m  ?TT\  \  b 

C(X)  =  £  Ck  cos  27rXk  A  =  X  ck e 


-m 


-m 


Then  the  convolution  of  C  and  W, 

I 

W  ‘ 


(X)  =  A  J  w  M  C  (X-cu)duj 
"  TE 


=  A 


??W 


_  27Ti  Xk  A  f  YK  27ri(j-K)a)A 


icke 


/ 


doi 


2zr 


Ickwk e 


27Ti  X  k  A 


is  the  modified  window.  This  equation  allows  W(X)  to  be  cal¬ 
culated  directly  in  terms  of  the  simple  window.  It  also  furnishes 
a  basis  for  selection  of  C(X)  or  ck  to  attain  a  desired  W{X) 

We  have  shown,  however,  that  we  do  not  have  complete  freedom 
in  selecting  a  window,  and  (2.  12)  can  only  be  a  guide  in  these  matters. 

The  modified  windows  in  usual  use  are  wider  and  have 
smaller  side  lobes  than  does  the  simple  window.  Graphs  of  some  of 
them  may  be  seen  in  Blackmand  and  Tukey,  or  in  Jenkins  (Dover,  and 
Technometrics,  v.  3,  2  referred  to  above).  In  general  they  compare 
to  an  ideal  rectange  centered  at  Xq  and  of  width  l/2m  A  .  if  this 
ideal  were  the  case  and  we  strung  out  such  rectangles  without  over¬ 
lapping  we  could  put  about  2m  +  1  in  the  interval  (-1/2A,  1/2A). 

Each  rectangular  window  would  "view"  a  separate  part  of  the  power 
spectrum  and  give  estimates  of  f (X)  at  2m  +  1  evenly  spaced 
points  in  the  interval  where  it  is  to  be  estimated.  If  we  take 
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the  spacing  of  the  points  and  the  width  of  the  rectangles  to  be 
l/(2m+l)  then  the  rectangles  are  somewhat  narrower  but  they 
do  not  go  outside  the  interval  where  f(X)  is  to  be  estimated. 

Let  us  use  the  notation  (see  (2.6)  and  (2.9)  ) 

m  m  l  j 

fK=  A  I  ricOS27T(kS)  jA  =  A  2  rj  COS  2  77-  — (2. 13) 

-  m  —  m  4rnt| 

for  the  spectral  estimates  at  evenly  spaced  points  using  the  sim¬ 
ple  window.  Then  by  using  either  the  exponential  form  of  this 
series  or  the  formula  cited  in  connection  with  (2.10) ,  we  have 
the  inverse  relation 

rj\SifKC0S27r7dr  <2.14) 

The  values  f^  are  sometimes  called  the  "raw  spectrum".  The 
modified  window  gives  the  spectral  estimates 

A  ™  ik 

fks  Z  rici  cos  27r  o  'Y  (2.15) 

K  _m  2m  +  l 

The  convolution  theorem  gives  (taking  f^  as  periodic  and  even) 

A  m 

fk  =  z  Vi  c  <i3)  (216) 

-m  ' 

where  C(j  6)  is  given  by  (2.11): 

m  ik 

C(jS)=  l  cKcOS27 

—  m 


(2.17) 
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With  evenly  spaced  spectral  estimates,  therefore,  the  modification 
of  the  window  may  be  done  by  a  running  average  smoothing  of  the 
raw  spectrum,  using  as  coefficients  the  finite  Fourier  transform 

of  the  c.|s. 

J 

The  inverse  relation  (2.13),  (2.14)  and  the  convolution 
relation  (2.  15),  (2. 16),  (2. 17)  are  to  be  used  later  in  a  discussion 
of  filtering. 

In  this  section  we  have  discussed  the  limitations  imposed 
by  the  finiteness  and  discreteness  of  the  data  on  what  we  must  try 
to  estimate.  We  shall  go  on  now  to  discuss  the  random  deviations 
from  this. 


3.  Statistical  problem 

We  have  been  considering  the  character  of  the  expected 
value  of  linear  spectral  estimates.  The  actual  outcome  of  the  es¬ 
timates  is  random,  the  ultimate  randomness  being  in  the  smoothed 
profile  data  itself.  We  will  assume  that  this  is  approximately 
Gaussian.  Evidence  in  support  of  this  assumption  will  be  found  in 
the  following  section.  We  discuss  here  the  main  ideas  of  the  sta¬ 
tistical  problem  without  giving  detailed  proofs.  These  may  be  found 
in  Blackman  and  Tukey  (Dover,  cited  above)  or  in  Grenander  and 
Rosenblatt  (Statistical  Analysis  of  Stationary  Time  Series,  Wiley 
1957).  In  Chapter  3  of  this  report,  the  section  corresponding  to 
this  one  has  similar  proofs  worked  out  for  the  two  dimensional  case. 

The  estimate  W  (2.6)  considered  as  a  function  of  the  n 
data  points,  is  a  homogeneous  quadratic  polynomial  (a  quadratic 
form).  In  practice  such  polynomials  in  Gaussian  random  variables 
are  found  to  have  distributions  that  are  well  approximated  by  a 
gamma  distribution,  i.e. 


W  <  2 


k  k  —  I 


roo 


-  At  . 
e  dt 


(3.1) 


-28- 


For  purposes  of  using  tabulated  values  it  is  convenient  to  recog- 

2 

nize  this  as  a  generalization  of  a  chi- squared  of  x  -distribution. 

v*-*  2 

In  fact  (3.  1)  would  be  stated  equivalently  as,  W/2A  has  a  x 
distribution  with  2k  degrees  of  freedom,  except  that  2k  need 
not  be  an  integer. 

The  gamma  distribution  is  fitted  by  the  method  of  moments, 
that  is  the  values  of  A  and  k  are  chosen  so  that  the  first  two  mo¬ 
ments  of  the  gamma  distribution  match  those  of  W. 

4  Ak  =  E  {w} 

16  A2k  =  Var  {w} 


from  which 
k 


A  = 


e{w} 

4k 


2 

Using  a  x  table  we  may  find  a  95%  confidence  interval  for 
E  {W}  by  getting  numbers  L  and  R  such  that 

P  {l<Xz<r}  =  .95 

Then  since  W/2A  =  2kW/E{W}  has  a  x^  distribution 


(3.2) 


p{-^-W<  E  {w}<  AA-w}  *  .95 


The  factors  2k /R  and  2k /L  give  the  ends  of  a  confidence  interval 
for  the  random  deviations  of  the  observed  W  from  the  expected  W. 
A  short  table  will  illustrate  how  the  width  of  the  interval  depends 
on  the  degrees  of  freedom,  2k. 
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2k 

2k  =  deg.  fdm.  R 

2k 

L 

5 

.  39 

6.  02 

10 

.  495 

3.  08 

15 

.  546 

2.  40 

20 

.  622 

2.  09 

25 

.  615 

1.  91 

30 

.  640 

1.79 

40 

.675 

1.  64 

50 

.  700 

1.  54 

60 

.  720 

1. 48 

The  degrees  of  freedom,  2k,  measure  the  amount  of 
statistical  error.  Various  equivalent  measures  of  statistical 
error  are  listed  in  the  paper  by  Jenkins  cited  above  (Technometrics 
v.  3,  No.  2).  They  all  depend  ultimately  on  the  first  two  moments 
of  W.  The  first  moment  of  W  is  given  in  equation  (2.7)  which 
we  repeat  now 

E  {w}  =  J  f  (X)  W(\)d\  (2.7  rep) 

-oo 

The  variance  of  W  may  be  approximated  in  a  similar  form  by 
making  use  of  the  assumption  of  approximate  normality. 

v°r  {") =  /ro,!<x>  w'lxidx 

-00 


(3.3) 
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The  proof  of  this  formula,  which  is  long  and  tedious,  may  be 
found  in  Blackman  and  Tukey  (Dover  1958,  cited  above)  and  will 
be  given  in  the  next  chapter  for  the  two  dimension  case.  With 
these  two  formulas,  the  degrees  of  freedom  formula  (3.  2)  becomes 


2k  =  2A  (n-m) 


[/f(X)  W(X)dX]2 
/  f2(X)W2(X)dX 


If  the  power  spectral  density  function,  f(A),  does  not  have  a  steep 
slope  in  the  region  where  the  window  is  large,  and  if  it  vanishes 
substantially  outside  the  interval  (-1/2  A>  1/2A)  then  the  integrals 
may  be  approximated  by  taking  an  average  value  of  f(X)  out  and 
these  average  values  will  approximately  cancel  giving 

[/¥  W(X)dX]2 

"2  a 

2  k  =  (  n  -  m  )  - ■ - 

/2,a  W2(X)dX 
“  IT 


The  numerator  of  this  should  be  unity  for  a  properly  normalized 
window  function,  so  the  size  of  the  denominator  gives  the  depend¬ 
ence  of  the  statistical  error  on  the  window.  Generally  it  is  larger 
for  the  sharper  windows.  For  an  ideal  window,  a  rectangle  of 
width  l/2m  A  and  height  mA  (to  make  a  unit  area  for  the  two  peaks; 
for  the  window  at  A  =  0  where  the  two  peaks  coincide,  the  height  should 
be  2m  A  .  )  we  have  a  final  approximation 


_l _ 

2  m2  A? 
2  A2 


n-m 


2k  -  2A(n  -  m) 


2 


m 
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This  approximation  is  on  the  optimistic  or  large  side  but  it  has  the 
advantage  of  simplicity. 

The  magnitude  of  the  statistical  error  may  be  indicated 
by  confidence  intervals.  These  may  be  found  from  the  degrees  of 
freedom  as  in  the  table  given  in  this  section.  From  the  table  one 
can  see  that  precision  increases  with  an  increase  in  degrees  of  free¬ 
dom.  Formula  (3.6)  suggests  that  degrees  of  freedom  increases 
as  the  ratio  g  =  m/n  decreases.  To  decrease  this  we  must  either 
raise  the  number,  n,  of  observations  or  lower  the  number,  m,  of  lags. 
The  former  almost  always  increases  the  expense  of  the  investigation, 
and  the  latter  decreases  the  resolution  of  the  spectral  estimates. 

4.  Smoothing  the  Profile 

In  the  first  section  we  immediately  passed  from  actual 
heights  h(x)  to  differences  d(x)  of  the  actual  heights  from  smoothed 
heights,  eqn.  (1.1).  We  postponed  discussion  of  this  until  the  present 
section  so  that  we  could  relate  this  discussion  to  the  ideas  of  the  in¬ 
tervening  sections.  Our  mathematical  assumption  concerning  the 
profile  of  ground  heights,  is  that  it  is  stationary- -of  the  same  char¬ 
acter  throughout  the  region  of  measurement.  In  most  areas  this 
property  is  clearly  untrue  for  any  possible  region,  though  it  might  be 
"maoe  true"  for  suitable  regions  if  the  hills  could  be  flattened.  In 
their  turn  the  hills  might  appear  as  roughness  in  very  large  regions 
in  which  general  statewide  level  changes  were  removed  and  so  on. 

Our  problem  is  to  carry  out  the  investigation  on  a  scale  relevant  to 
our  purpose. 

The  first  part  of  limiting  the  scale  of  the  investigation  is  to 
limit  the  extent  of  the  region  to  be  measured  so  that  the  roughncss-- 
other  „han  "hils"  is  of  a  uniform  character.  The  second  part  is  to 
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remove  the  hills  in  some  way.  In  the  limited  region  they  are  not 
usually  stationary.  (In  the  Battlefield  Day  data  the  ground  rises 
steadily  so  that  the  heights  at  the  end  are  not  nearly  the  same  as 
those  at  the  beginning. )  They  are  also  not  relevant  as  vibration 
producing  roughness.  Nevertheless  the  variations  in  height  due  to 
these  "hills"  are  ordinarily  much  larger  than  those  variations  we 
want  to  study.  If  we  were  to  leave  them  in  the  data  the  spectrum 
would  consist  almost  exclusively  of  a  central  value  at  zero  fre¬ 
quency.  The  other  values  would  be  lost  in  the  round-off  of  the  com¬ 
putation. 

We  have  used  two  methods  of  smoothing  the  profile.  One 
method  consists  in  constructing  a  smooth  profile  by  fitting  para¬ 
boloids  to  segments  of  the  data.  The  second--so  far  in  use  only 
with  the  linear  data--is  a  running  average  scheme. 

The  first  method  consists  in  building  up  a  smoothed  pro¬ 
file  by  fitting  second  degree  polynomials  F(x)  to  the  data  points 
h  (x^).  Actually,  several  polynomials  may  be  fitted  each  to  a  sep¬ 
arate  segment.  We  may  imagine  these  to  be  graphed  along  with 
the  data  so  that  the  data  or  actual  ground  profile  cuts  back  and  forth 
across  the  smooth  parabolical  curve.  It  is  clear  that  a  parabola  can 
remove  only  one  hill  from  each  segment  so  that  in  a  profile  with 
several  hills  there  should  be  several  segments  fitted  by  parabolas. 
This  was  done,  but  it  was  found  that  the  smoothed  profiles  obtained 
in  this  way  hau  steps  where  the  segments  came  together.  In  some 
of  the  data  the  steps,  though  few  in  number,  were  about  ten  times  the 
values  of  the  d(x)  values  and  an  approximate  analysis  indicated 
that  they  could  affect  the  resulting  spectra.  (See  Appendix  Y  of 
this  chapter. 
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The  theoretical  effect  of  parabolic  smoothing  is  difficult 
to  work  out.  Computations  were  made  which  indicated  that  this 
type  of  smoothing  has  several  advantages  over  the  running  average 
type  of  smoothing.  We  have  experienced  some  difficulties  in  the 
analysis  of  our  data  which  suggest  that  parabolic  smoothing  works 
well  when  the  parabolas  fit  closely  but  stronger  smoothing  is  needed 
on  hilly  ground. 

We  have,  therefore,  done  some  calculations  using  another 
spectral  program  which  has  been  equipped  with  a  smoothing  sub¬ 
routine  of  the  running  average  type.  A  theoretical  discussion  of  the 
effects  of  this  type  of  smoothing  is,  fortunately,  very  easy  and  we 
will  present  it  for  the  linear  case. 

Let  us  denote  the  heights  of  the  ground  above  points  evenly 
spaced  on  a  straight  line  of  the  datum  plane  by 


-a  ♦  l 


i  •  •  •  i  h.  .  h. 


^*n  ^  n  ♦!  ’  * 


n*  a 


(4.  1) 


We  have  added  points  on  each  end  for  a  purpose  which  will  be  clear 
later.  Running  average  smoothing  requires  that  we  select  some 
coefficients 


*  k-a  ♦!••••»  b©  *  •••»  bg  (4.2) 

The  smoothed  values  will  be  calculated  by  the  formula 

o 

8  £  hk*a  k  =  1 , 2  , .  •  • »  n  (4.  3) 

a=-a 
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(Because  of  our  foresight  the  numbers  here  run  from  1  to  n.  Sim¬ 
ilar  foresight  in  the  computations,  accomplished  by  faking  in  a  few 
values  at  the  beginning  and  end  of  the  data  accomplishes  the  same 
end. )  The  deviations  used  for  spectral  calculation  and  called  "the 
smoothed  profile"  are: 

o  a 

bk  ^k  "  bk  2  ba  hk+a  "  £  ba  bk*a  (4.4) 
Cf=-a  a=_a 

where 


-ba 
b’  =  a 

if 

a  =  o 

a 

<‘bo 

for 

a  =  o 

To  show  what  this  type  of  smoothing  can  accomplish,  let 
us  suppose  that  the  ground  heights  (5.  1)  are  composed  of  a  smooth 
"hilly"  part  plus  a  rought  part  6^ 

hk  *  Pk  +  ^k 

Then 


dk  "  £  ba  ^k*a  +  £  ba  ^k«-a 

Now  if  the  P^  are  values  of  a  polynomial  on  integer  points 

Pk  =  C0+C,  k  +C2  k2  +  •  •  •  Cmkm 


(4.5) 
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Then  the  coefficients  b  (4.  2)  may  be  selected  so  that  the  first 
term  of  (4.5)  is  zero.  Suppose,  for  example,  the  is  of  2nd 
degree .  Then 


£  ba  Pk*a 
a  =-a 


a=-a 


0+C,(k+a)  +  c2(k+ 


■  S  b’  |(C0+C,  k+C2  k2)  +  a  (C,  +  2C2k)  +  a2  C2  I  (4.6) 
a  L  J 


=  (C0+C,  k+  C2k2)  I  b*  +  (C,+2C2k)  lab’  +  C2  Za2b’ 

a  a  a 


If  we  select  the  coefficients  so  that 


I  ab^  =0 

a 


aVa=0 


(4.7) 


Then  the  right  side  of  (4.  6)  will  be  zero.  In  this  sense  the  running 
average  smoothing  can  remove  a  parabolic  trend.  Note  that  the 
smoothing  is  very  strong  in  that  the  polynomial  P  may  be  different 
at  each  point.  Since  (4.7)  is  a  set  of  three  homogenous  linear  e- 
quations  with  coefficients 

1,  1,  .  ..  1,  1,  1,  •  •  •  1,  1 
-a,  -a+1,  .,-1,  0,1,..  .  a-1,  a 
a2,  (a-1)2,  .  .,1,0,1,  ..(a-1)2,  a2. 


For  a  =  0  or  1  the  matrix  is  non  singular  so  there  is  only  the 
trivial--all  zero- -solution.  Thus  a  must  be  at  least  2.  The  so¬ 
lutions  are  not  unique  and  for  larger  a's  there  is  a  large  variety 
of  them. 

The  removal  of  trend  is  not  the  only  effect  of  this  smooth¬ 
ing.  In  (4.  5)  the  second  term  on  the  right  is  the  smoothed  time 
series  whose  spectrum  we  want.  When  we  estimate  the  spectrum  of 
the  smoothed  series  we  will  not  be  estimating  what  we  want.  Sup¬ 
pose  we  denote 

=_2  ba  Sk  +  a  k  =  1,2 . n 

The  covariance  function  for  7^  is 

rd  s  E  {  ^k  ♦  d  }  =  E  {  I  ba8k+a  I 

Ct=-a  /?  =  -a 

=  |  |  babi  E(8k.a8h.d.„) 

=  a  is  ^  W*a 


where  r  =  E{6,  6, +a  }  is  the  covariance  function  of  the  6's.  If 

cl  K.  K 

we  change  variables  in  the  summation  so  that  (3  -  a  =  p  and  con¬ 
sider  the  b's  to  be  defined  for  all  a,  taking  the  value  zero  for 
a  <-  a,  a>  a  then  the  covariance  function  of  the  smoothed  series 
may  be  written 
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which  is  a  running  average  smoothing  of  the  y's  using  the  coefficients 


(4.9) 


A  running  average  smoothing  of  the  stationary  part  of  a 
profile  results  in  a  profile  whose  covariances  are  a  running  average 
smoothing  of  the  original  covariances.  The  coefficients  are  the 
"covariances"  of  the  coefficients  used  to  smooth  the  profile. 

The  running  average  (4.8)  is  a  convolution.  If  we  rewrite 
the  convolution  relations  (2.15),  (2.16),  (2 .17)  using  the  same  con¬ 
ventions  as  to  periodicity,  we  have 
~  m 

rd  =  Z  rd+j  Bj  (2. 16  rep) 

j  =  -m 


where 


jk 

2  m  +1 


m 

/3:  =  Z  Bk  COS  27T 

1  —  m 


jk 

2  m+l 


(2 . 15  rep) 


and  f.  is  the  r's  given  by  (2.14). 

3 

A  running  average  smoothing  of  the  stationary  part  of  a 
profile  results  in  a  profile  whose  spectral  estimates  are  multiples 
of  those  of  the  original  profile.  The  factors  are  the  "power  spec¬ 
tral  estimates"  of  the  smoothing  coefficients.  This  is  so  since 

(2.15  rep)  shows  the  spectrum  of  r  ,  to  be  f.j3..  If  we  calculate  f., 

^  J  J  J 

the  spectral  estimates  of  the  smoothed  data,  then  we  may  correct 
it  using  as  factors  1  /  fi.,  given  by  (2.17  rep). 
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Some  work  remains  to  be  done  on  this  type  of  smoothing 
that  needs  experience  more  than  mathematical  analysis.  The 
smoothing  coefficients,  b',  which  we  have  tried,  have  a  spectrum 
which  is  zero  at  X  =  0,  and  very  close  to  zero  for  low  frequencies. 
The  extent  of  the  frequency  range  for  which  these  spectral  values 
are  zero  indicates  the  severity  of  smoothing  by  indicating  how  much 
of  the  stationary  part  of  the  profile  is  removed.  Correcting  may 
not  give  it  back  since  it  also  exaggerates  round-off  error  as  well  as 
the  low  spectral  value.  Therefore,  some  care  must  be  exercised 
in  choosing  the  b'  so  as  to  leave  in  the  important  part  of  the  desired 
spectrum.  We  may  choose  from  the  variety  of  smoothing  coefficients, 
available  to  us  as  solutions  of  equations  such  as  (5.7),  those  which 
give  a  desirable  spectrum --equivalently  a  desirable  covariance  func¬ 
tion  (5.  8).  In  general  there  may  be  many  solutions  of  (5.  7)  giving 
exactly  the  same  covariances.  Each  of  these  will  act  the  same  ac¬ 
cording  to  our  mathematical  model;  they  will  remove  the  appropriate 
polynomial  trend  and  affect  the  spectrum  of  the  stationary  part  in 
the  prescribed  way.  What  is  the  actual  effect  in  practice  depends  on 
how  our  actual  problem  differs  from  this  mathematical  problem. 

This  will  require  an  investigation  by  experience  rather  than  by  logic. 

The  problem  of  finding  coefficients  that  have  a  prescribed 
spectrum  i.e.  ,  prescribed  correction  factors  seems  to  be  easier  to 
solve  by  computation  than  by  analysis.  This  is  another  area  for  ex¬ 
perience.  As  yet,  we  have  used  only  one  set  of  coefficients  so  these 
matters  are  for  future  investigation. 


Inasmuch  as  the  coefficients 


b'. 

J 


may  be  chosen  so  as  to 


make  the  spectrum  f^  have  smoother  slopes,  it  will  help  in  the 


estimation  problem.  There  will  be  no  sharp  bumps  that  leak  in 
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through  the  side  bands  and  so  on.  Smoothing  for  this  purpose 
has  been  called  "prewhitening".  Our  first  computer  program  has 
an  elaborate  prewhitening  section  in  which  the  spectrum  was 
first  estimated  and  then  coefficients  were  calculated  to  correct 
it.  Our  experience  indicates  that  the  main  bad  part  of  the  spec¬ 
trum  occurs  at  zero  frequency  so  that  simpler  methods  of  pre¬ 
whitening  are  more  effective. 

5.  The  Computing  Programs 

We  have  used  two  computer  programs  to  estimate  line 
spectra.  A  complete  description  of  the  first  program  is  no  longer 
available.  The  general  operation  of  that  program  is  as  follows: 

1)  The  data  is  "smoothed"  by  taking  deviations  of  the  pro¬ 
file  heights  from  heights  of  quadratic  polynomials  fitted  to  segments 
by  least  squares. 

2)  Covariances  (mean  lagged  products)  are  computed. 

3)  A  preliminary  spectrum  is  calculated  by  making  a 
cosine  transform  of  the  covariances  and  smoothing  it  by  the  Ham¬ 
ming  coefficients. 

4)  A  trigonometric  polynomial  of  a  given  number  of  terms 
is  fitted  to  the  spectrum  by  least  squares.  The  coefficients  of  this 
are  used  in  a  series  of  equations  whose  solutions  are  running  average 
coefficients  to  be  used  on  the  profile  data.  This  running  average 
affects  the  spectrum  as  outline  din  Section  4  and  the  computations  are 
made  so  that  it  tends  to  flatten  out  steep  places  of  the  spectrum. 

This  is  the  "Prewhitening"  referred  to  in  Section  4. 

5)  The  profile  is  smoothed  by  running  average. 

6)  New  covariances  are  computed  and  also  their  spectrum 
(cosine  transform). 
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7)  Correction  coefficients  for  the  prewhitening  running 
average  are  computed. 

8)  A  corrected  spectrum  is  computed. 

9)  The  corrected  spectrum  is  smoothed  by  the  Hamming 
coefficients. 

This  program  recently  became  inoperable  and  we  have  re¬ 
placed  it.  We  have  obtained  a  spectral  program  made  for  the  IBM 
7090.  The  complete  details  of  the  program  are  available.  Generally 
it  performs  these  operations: 

1)  Subtracts  the  average  of  the  entire  profile  from  each 
profile  value. 

2)  Computes  covariances  of  any  desired  segment  of  the 
profile  data. 

3)  Computes  the  spectrum  of  the  covariances  using  the 
Hamming  W  function. 

The  first  step  is  a  sort  of  smoothing  but  it  is  not  strong 
enough  for  our  purposes.  We  have  modified  the  program  by  insert¬ 
ing  a  subroutine  which,  previous  to  the  above  steps 

(0)  smooths  the  profile  by  running  average  using  specified 
coefficients. 

This  modification  is  only  temporary  and  it  is  lacking  in  several  re¬ 
spects.  The  most  serious  respect  is  that  the  correction  factors 
must  be  calculated  and  applied  by  hand.  Since  this  type  of  smooth¬ 
ing  seems  to  be  very  satisfactory,  we  are  including  a  new  sub¬ 
routine  which  will 

(0.  0)  smooth  the  profile  by  running  average  using  specified 
coefficients. 

(0.  1)  indicate  any  profile  values  that  are  inconsistent  with 
their  neighbors. 

(0.2)  compute  and  store  the  correction  factors 

(0.3)  correct  the  spectral  estimates. 
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6.  Discussion  of  One -dimensional  Spectral  Results 

In  this  section,  we  shall  comment  in  some  detail  upon 
the  one -dimensional  p.  s.  d.  's,  how  they  were  obtained,  and  their 
corresponding  profiles. 

One  of  the  significant  problems  from  the  point  of  view  of  sub¬ 
sequent  use  is  the  taking  out  or  elimination  of  trends,  hills,  etc. 
which  are  not  considered  relevant.  When  we  began  the  survey  data 
analysis,  we  estimated  trends  by  means  of  quadratic  polynomials 
fit  to  segments  of  the  data  by  least  squares.  These  polynomials 
were  subtracted  from  the  profile  heights  to  produce  smoothed  data 
(deviations  from  fitted  polynomials).  The  roughness  of  these  smooth¬ 
ed  data  were  then  analyzed  in  the  spectral  analysis  program.  The 
spectral  estimated  for  the  Yuma  data  contained  anomalies  such  as 
high  tail  values  and  occasional  negative  values.  In  the  case  of 
Aberdeen  and  Knox,  these  anomalies  were  not  observed.  Moreover, 
in  the  Yuma  data,  spectral  estimates  of  the  halves  and  whole  were 
different  in  character.  It  was  not  possible  to  interpret  these  anom¬ 
alies  in  terms  of  the  data.  This  was  somewhat  surprising  in  view  of 
the  consistancy  of  the  results  with  Aberdeen  and  Knox  and  from  high¬ 
way  survey  data. 

The  origin  of  these  anomalies  remained  a  mystery  for  some 
time,  but  attention  was  gradually  directed  to  the  smoothing  proceedure. 
We  next  obtained  another  p.  s.  d.  estimate  program  which  only  re¬ 
moved  the  general  average;  it  produced  with  the  smoothed  Yuma  data 
precisely  the  same  p.  s.  d.  estimates  which  established  beyond  doubt 
the  troubles  in  the  smoothing  proceedure.  We  then  made  a  prelim¬ 
inary  routine  for  running  average  smoothing  which  added  to  the  pro¬ 
gram  produced  p.  s.  d.  estimates  free  of  most  of  the  anomalies  noted 
above.  Moreover,  the  remaining  anomalies  were  immediately  trace¬ 
able  to  inconsistant  raw  data  points  not  previously  detected.  The 
results  presented  here  are  those  obtained  from  this  revised  program. 
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The  general  features  were  already  discussed  in  the  Intro¬ 
duction  and  that  discussion  will  it  be  repeated.  We  must  describe 
however  the  significance  of  data  processing  and  analysis. 

As  explained  in  Section  II,  Part  4,  smoothing  throws 
away  part  of  the  data  and  partially  restores  it  by  numerical  cor¬ 
rection.  The  nature  and  significance  in  vehicle  dynamics  of  the 
part  thrown  away  is  not  completely  known  to  us  and  the  control  of 
this  part  by  choice  of  smoothing  coefficients  is  not  yet  entirely  com¬ 
plete.  It  is  important  to  understand  that  smoothing  must  bare  a 
direct  relationship  to  p.s.d.  estimate  use.  Whether  hills  are  re¬ 
moved  or  not  depends  upon  vehicle  size,  speed  and  the  presence 
and  nature  of  the  hills.  Thus,  use  of  the  p.  s.d.  estimate  procedures 
and  the  estimates  without  regard  to  the  method  of  estimation  em¬ 
ployed  may  in  unusual  conditions  lead  to  error. 

Figure  5  shows  the  p.  s.  d.  estimates  of  the  full  750  data 
points,  the  first  350  data  points,  and  the  second  350  data  points  from 
Aberdeen.  The  raw  data,  the  smoothed  data,  the  covariance  of  the 
smoothed  data,  the  p.  s.d.  estimates  of  the  smoothed  data,  and 
the  corrected  p.s.d.  estimates  are  given  in  Tables  1  to  ID.  The 
plots  are  of  the  corrected  p.  s.  d.  estimates.  The  uncorrccted  p.  s.  d. 
estimates,  calculated  from  the  smoothed  data,  arc  distorted  by  the 
smoothing  and  are  corrected  by  multiplication  as  explained 
in  Section  II,  Part  4. 

In  general,  the  p.  s.  d.  of  the  full  set  of  data  should  be  ap¬ 
proximately  the  average  of  the  p.  s.  d.  's  of  the  two  halves.  This 
feature  is  noticeable  upon  Figure  5.  The  matching  in  the  significant 
range  \  =  .005,  to  X  =  0.075  between  the  halves  is  substan¬ 
tial  except  in  the  interval  containing  the  bump.  The  bump  is  due 
to  a  periodic  component  in  the  profile,  (Figure  1)  which  was  most 
visible  in  the  first  half  of  the  profile.  This  is  confirmed  by  the  pres¬ 
ence  of  a  strong  bump  in  the  first  half  p.  s.  d.  and  the  absence  of  any 
bump  in  the  second  half.  The  generally  high  level  in  the  low  tails 
is  due  to  inconsistencies  in  the  data  as  evidenced  by  the  readings  at 


-43- 


112  feet  and  1,254  feet  in  the  smoothed  data  (Table  IB).  As  these 
tails  are  not  of  interest,  these  inconsistencies  are  not  important 
to  us  at  this  time.  The  other  significant  difference  between  halves 
is  in  the  X  =  .09  toX  =  .10,  where  the  second  half  seems  to  have 
about  twice  as  much  power  as  the  first.  Examination  of  the  profile 
to  confirm  this  by  inspection  indicates  that  visual  methods  are  not 
capable  of  explaining  this  difference,  but  since  the  power  involved 
is  low  this  is  not  surprising. 

Figure  6  gives  the  p.  s.d.  estimates  of  the  full  and  two 
half  sets  of  data  points  from  Ft.  Knox.  The  tabulated  results  are 
given  in  Tables  2  to  2D. 

The  three  spectra  match  better,  in  general,  than  those  of 
Aberdeen.  In  particular,  the  tails  fall  off  far  more  rapidly  and 
to  a  lower  value  than  in  Aberdeen.  The  broad  bump  in  the  first 
half  p.  s.  d.  is  spread  over  the  wave  lengths  from  25  feet  to  50  feet 
whereas  in  the  second  half  this  power  is  concentrated  more  in  the 
wave  length  around  30  feet.  Again  these  differences  are  not  notice¬ 
able  in  scanning  the  profile  (Figure  2)  by  eye. 

Figure  7  gives  the  3  p.  s.d.  estimates  from  the  Yuma 
data.  The  tabulated  results  are  given  in  Tables  3  to  3D. 

The  matching  in  Yuma  is  quite  remarkable.  If  anything, 
the  first  half  p.s.d.  estimate  indicates  less  roughness  in  the  higher 
frequencies  than  the  second  half;  this  is  evident  in  the  profile 
(Figure  3).  The  drop  off  to  a  very  low  level  of  power  in  the  tails 
is  also  very  satisfactory. 

Figure  8  gives  the  5  p.s.d.  estimates  from  the  Battlefield 
Day  data.  Numerical  data  and  computed  results  are  tabulated  in 
Tables  4  to  4D.  In  this  figure,  we  give  the  p.  s.d.  estimates  of  the 
four  quarters  as  well  as  of  the  whole.  The  estimates  match  verv 
well  up  to  X=  .08.  There  is  no  noticeable  bump  or  any  other  feature 
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of  interest  in  the  range  from  X.  =  0.00  to  X  =  .  08.  Two  of  the  quar¬ 
ters  have  comparatively  high  tails  and  one  exhibits  a  distinct  perio¬ 
dicity.  These  features  are  transmitted  to  the  spectra  estimate  of 
the  full  set  of  data.  They  are  adequately  explained  by  the  presence 
of  inconsistent  data  points,  as  mentioned  in  Appendix  Y,  Section  2. 

The  differences  between  p.  s.  d.  estimates  of  the  different 
tracts  of  ground  are  much  greater  than  the  differences  encountered 
with  estimates  from  sections  of  the  same  tract.  This  indicates  that 
the  assumption  of  stationarity  over  reasonable  lengths  is  a  good  one 
and  justifies  the  use  of  p.  s.d.  techniques  in  the  measurement  of 
ground  roughness.  It  must  be  expected  that  the  assumption  of  sta¬ 
tionarity  is  not  ideally  fulfilled  in  any  tract  of  ground.  However, 
in  ground  selected  with  reasonable  care,  our  data  reveals  that  the 
assumption  is  reasonably  well  fulfilled  over  lengths  of  1,000  feet 
or  more.  This  indicates  the  possibility  of  adjusting  suspension  for 
slowly  varying  roughness. 

These  p.s.d.  estimates  have  similar  shapes  though  they 
are  from  widely  different  areas  of  the  country.  The  fact  that  the 
spectral  shapes  are  similar  and  simple  in  form  is  a  matter  of  some 
importance.  First,  it  might  be  possible  to  character  ground  rough¬ 
ness  (spectral  shape)  in  terms  of  three  or  four  parameters.  This 
may  simplify  classification  of  ground  roughness  types;  it  would  per¬ 
mit  easy  simulation  of  roughness  types;  and  it  would  define  ranges 
of  suspension  system  parameters  for  adjustment  to  changing  ground 
roughness. 

Figures  9,  10,  and  11  show  cumulative  diagrams  of  the 
smoothed  data  from  Aberdeen,  Knox,  and  Yuma  plotted  on  proba¬ 
bility  paper.  The  scale  is  the  same  in  each;  the  horizontal  scale 
is  in  percentage  of  the  data,  the  vertical  in  feet  of  deviation.  The 
points  are  plotted  to  indicate  the  percentages  of  data  smaller  than 
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various  values.  The  arrangement  of  values  on  the  horizontal  scale 
is  such  that  a  theoretical  Gaussian  distribution  plots  as  a  straight 
line  passing  through  the  points 

(  15.87,  jJL- a)  ,  (0  ,  fJL)  ,  (84. 13, ft  +  C7  ) 

2 

where  n  is  its  mean  and  cr  its  variance.  Gaussian  data  will  also 
be  straight  except  for  random  deviations.  The  values  of/i  and  cr 
may  be  estimated  by  fitting  a  straight  line  and  reading  the  ordinates 
of  the  points  15.87,  0,  84.13.  Interpretation  of  deviations  from 
straightness  requires  some  experience  but  is  not  difficult. 

For  a  random  function  satisfying  the  ergodic  hypothesis, 
the  time  average  of  an  indicator  function,  I  (A),  which  is  unity  if 
the  event  A  occurs  and  zero  if  not,  equals  the  ensemble  average 
or  probability  of  A.  This  means  for  profile  heights  that  the  fraction 
of  heights  smaller  than  a  number,  z,  should  approximate  the  prob¬ 
ability  that  a  height  is  smaller  than  z.  If  the  profile  is  a  Gaussian 
random  function,  the  data  should  therefore  be  Gaussian  data.  The 
curves  for  our  data  indicate  systematic  departures  from  straight¬ 
ness.  In  our  opinion  the  amount  of  departure  is  not  enough  to  in¬ 
validate  the  Gaussian  approximations  we  must  use. 

The  curves  suggest  that  the  distribution  of  the  data  is 
symmetric,  as  is  the  Gaussian  distribution,  but  that  it  has  a  higher 
percentage  of  its  probability  in  the  region  of  small  deviations  than 
does  the  Gaussian.  This  might  possibly  be  a  result  of  the  effects  of 
erosion  to  round  off  peaks  and  fill  in  holes.  The  amount  of  this  de¬ 
parture  from  the  Gaussian  is  more  for  Yuma  than  for  Aberdeen  or 
Knox.  It  is  possible  that  we  may  find  in  the  future  that  the  deviation 
from  Gaussian  for  some  locations  is  large  enough  that  we  should 
inquire  into  its  effect  on  the  Gaussian  approximations.  This  effect 
will  be  confined  to  the  formulas  concerning  the  statistical  errors. 
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Appendix  Y 

Effect  of  Certain  Errors  on  Spectral  Calculations. 

In  this  appendix  we  will  present  the  results  of  two  inves¬ 
tigations  we  carried  out  in  an  attempt  to  explain  some  anomalous 
features  of  the  spectral  calculations  for  Yuma.  The  first  of  these 
was  done  after  one  of  the  profile  heights  was  found  to  be  some 
thirty  feet  lower  than  its  neighbors.  It  was  extended  somewhat  in 
the  light  of  similar  errors  found  in  other  profiles.  The  second  was 
concerned  with  the  effect  of  a  bad  joint  between  parts  of  a  profile 
that  has  been  smoothed  by  fitting  parabolas  to  two  segments.  The 
result  of  the  latter  indicated  that  the  unruly  behavior  of  the  Yuma 
spectrum  could  not  be  due  principally  to  this  type  of  error  so  the 
investigation  was  not  developed  further. 

For  the  first  investigation  we  will  suppose  that  the  proper 

ground  profile  data  djdg.-.d  is  actually  presented  to  the  computer 

as  d'.d'  d'  where 
1  2...,  n 


d’. 


=  dk 

=  da  +  A 


k  a 
k  =  a 


(Yl) 


The  covariances  actually  computed  are 

I  V'  _i  *  j  • 


=  TT^T  £  dkdk+t 

=  TTTf  I  dkdk+t  +  (  da+t  +  do-i  ) 

k 


(Y2) 


for  t  ^  0  ,  and 
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For  large  A  the  difference  between  V  and  the  neighboring 

2  0 

7*1,  71  g  etc.  (on  the  order  of  A  /n)  is  noticeable.  If  a  is 
near  one  end  of  the  data,  the  point  where  d  +  .  cuts  out  of  Y  2 
may  be  noticeable  (to  avoid  complicated  limits  of  summation  we 
take  =  0  for  k<0  or  k>n.) 

The  computed  spectral  values  are 


fu  + 


2  Ada 
n 


m 


T=  I 


dp+t  +  dp-t 


cos  2  7 r 


_ut _ 

2  m  ♦  I 


(Y4) 


In  this  formula  the  terms 


dp+t  +  d0.t 
n-t 


cos  2  7 r 


ut 

2mtl 


are  a  Fourier  series  of  a  small  segment  of  the  data  and  may  be 

2 

quite  erratic.  The  dominant  term  in  (Y.4)  is  A  /n  which  is  a 
constant  or  "white  noise"  term.  This  was  a  very  prominent  feature 
of  the  spectrum  obtained  from  the  Yuma  data  when  it  contained  a 
single  large  error  in  the  profile. 

A  single  small  error  will  have  its  size  divided  by  n,  the 
length  of  the  profile  record,  in  the  Fourier  series  and  its  square 
divided  by  n  in  the  white  noise  term.  Since  n  is  usually  large, 
a  small  error  will  not  have  much  effect.  The  effect  of  several  er¬ 
rors  may,  however,  accumulate  to  add  noticeable  white  noise  to 
the  calculated  spectrum. 


-48- 


The  third  quarter  of  the  Battlefield  Day  profile  data  pre¬ 
sented  a  spectrum  that  was  affected,  seemingly,  by  both  a  white 
noise  and  a  periodic  component.  The  latter  had  a  period  of  almost 
exactly  six  stations,  or  about  .06  cycles  per  foot.  Since  1/.06  = 
16.6  feet  one  can  expect  a  pip  on  the  covariances  near  station  17. 
Such  a  pip  is  indeed  to  be  found  there  as  well  as  a  pip  at  zero  to 
account  for  the  white  noise.  The  full  run  spectrum  also  had  a  sim¬ 
ilar  shape  and  has  similar  pips  so  these  are  not  machine  errors. 

The  data  was  searched  for  bad  points  to  explainthe  pip  at 
the  zero  covariances  and  the  white  noise  part  of  the  spectrum.  A 
list  of  those  found  is 


No.  317  5 
4474 
4540 
4557 


2.0  ft  above  adjacent  points 
5.0  ft  above  adjacent  points 
3.0  ft  above  adjacent  points 
10.0  ft  above  adjacent  points 


The  last  two  are  exactly  17  stations  apart  and  are  without  doubt 

the  cause  of  the  pip  at  the  17th  lag.  The  white  noise  parts  of  the 

errors  add  together,  as  we  shall  show,  and  adding  up  the  values 
2 

A  In  for  the  above  we  have 


(  4  +  25  +  9  +  lOO)  /  1550  =  .09 


The  height  of  the  zero'th  lag  for  the  third  quarter  data  is  .0916, 
white  the  first  lag  is  -.0112  making  a  pip  of  .0928  which  is  quite 
close  to  the  above  estimate. 

If  two  errors  are  farther  apart  than  the  maximum  lag  then 
their  effects  in  (Y2)  and  (Y3)  are  additive.  However,  if  two  errors,  say 
of  A,  and  A9  are  separated  by  a  distance  of  t  <  m  stations  in  the 
data ,  say 


*0  =  Vo  +  A. 


+  A 


2 
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Then  the  effects  of  these  errors  on  the  covariances  are 

’  r*  +  n  - 1 1 1  y<J+t +  ya-t  ^  +  ^2  l  y0+Vt  +  yo*t0-t^  } 
W  HTTV  W  W  V.0>  +A2(y0,vl+  y„.v,)}+ 


A^+A22 


r0  =  r0  +  n  (2A|ya  +  2A2ya  .  )  +  _! - 


The  effects  are  additive  as  for  widely  separated  errors  except  at 
a  lag  of  t  where  there  is  a  pip  of 

A|  a2 

n-  |  t0| 


For  the  Battlefield  Day  data  third  quarter  errors,  this  is 
3x10  /  (1550  -  17)  =  .0195 

In  the  covariances  the  17th  is  .0178  while  the  neighbors  are  -.0020 
making  .0198.  These  anomalies  in  the  spectrum  are  therefore  ex¬ 
plained  adequately  by  the  errors  or  inconsistencies  in  the  data. 

The  spectra  for  Ft.  Knox  and  Yuma  seem  to  be  affected 
somewhat  by  periodic  components  and  the  spectra  for  Aberdeen 
seems  to  be  rather  high  on  the  high  frequency  end.  The  compilation 
of  this  data  for  the  purpose  of  plotting  on  probability  paper  shows 
several  outlying  in  each  profile,  but  none  as  spectacular  as  those 
found  in  the  Battlefield  Day  data.  The  cumulative  effect  of  these 
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could  explain  the  behavior  of  the  spectra.  In  these  cases  the  actual 
amount  of  the  bad  behavior  is  small  and  seems  to  affect  mostly 
those  frequencies  that  are  unimportant  in  our  work. 

The  second  investigation  is  concerned  with  a  smoothed 
profile  d^dg..^  ^n+i“*^2n  are  deviations  of  the  actual  pro¬ 

file  heights  from  parabolas  fitted  to  each  half.  We  will  suppose  that 
the  fitted  parabolas  do  not  meet  in  the  middle--thus  the  d's  also 
do  not  meet--and  we  will  try  to  estimate  the  effect  of  this  error. 
Our  attention  will  be  on  the  relation  between  the  spectrum  f  '(u)  of 
the  first  half,  f  "  (u)  of  the  second,  and  f  (u)  of  the  entire  profile. 
In  theory  the  average  of  the  two  parts  should  give  the  whole.  To 
show  this,  let  us  recall 

I  2n-r 

rt  =  2n-t  £  dk  dk*r  f  >  0 

k=l 


We  will  take  t^O  for  the  sake  of  having  definite  limits  on  the  sum. 
This  sum  may  be  partitioned 


r.  = 


t  '  2n- 


r  13.  Zn-t  , 

V  [  Zdkdk*t+Z  dkdk+t  +  Z  dk  dk+1 1 


k=n-UI 


2n-t 

Z 

k=n+l 


2  n-t 


[(n-t)  ri 


k=n-t»l 


dk  dk*t  + 


(  n-t )r”  ] 


(Y5) 


r’  +  r” 
rt  rt 


I 


2  n-t 


I.Vk., 


n-t= 
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If  t  is  a  small  fraction  of  n-t  and  2n-t  is  large  with  respect  to  the 
sizes  of  the  data  near  the  joint,  then 


.  +  r" 

f  "  2 


(Y6) 


Since  the  cosine  transform  is  a  linear  operation  this  equation  leads 
to 


fu  =  \  (  fu  +  fu) 


The  approximation  to  the  average  in  the  first  term  on  the 
right  of  (Y.5)  may  be  made  somewhat  better  by  choosing  a  value  of 
t  near,  say,  l/3m  since  the  r’s  are  largest  near  zero.  With  this 
(Y.7)  may  be  modified  a  little 

=  (•--£-)  t  k+o 

The  central  values,  however,  are  not  as  predictable.  The  sum  is 
the  same  kind  as  used  to  estimate  pt  and  if  it  weren't  for  the  error 
we  are  discussing  we  might  expect  the  second  term  to  have  the  ap¬ 
proximate  size 

t  r  ~  I  r 

TT^T  rt  ~  TF  rt 

which  would  cancel  out  the  l/6n  part  of  (Y.8)  on  transformation. 

(Of  course  this  is  approximate  and  the  amounts  1/Gn  are  very  small 
--no  doubt  the  other  approximations  make  larger  effects  than  this.) 
The  error  in  matching  the  two  halves  of  the  profile  affects  this  central 
part  only. 
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We  will  assume  that  the  values  are  roughly  constant 
on  each  side  of  the  half  mark  and  are  separated  by  a  distance  D 


d  n  -  m  <  k  <  n 

*■  d  +  D  n  +  I  <  k  <  n  +  m 


where  m  is  the  maximum  lag.  Then  for  any  value  of  t  0  t  <Jm 
the  central  term  is 


_J _ 

2  n-t 


X 


k=n-t+l 


dk  dk  +  t 


t 

2  n+t 


d  (  d  +  D) 


d  (d  +  D)  t 

2  n 


Since  rt  is  an  even  function  we  have  for  -m^t^m 


rt  +  r"  .  d  {  d  +  D) 

f  -  - 5*  -  + 


2  n 


1 1 


The  cosine  transform  of  |t|  maybe  gotten  in  closed  form  in  formula 
(428)  of  Jolley's  "Summation  of  Series"  (Dover,  1961) 

(  m  +  t)  sin  IT  u  '  2 m  +  1 


X  I  t  I  COS  2  TT-^-- 

t=  -m 


.  m  ♦  I 

|  -  cos  — - -  2  7T  U 


sin 


"TT  U 
2m  ♦  I 


2  sin' 


7 Tu 
2  m  +  l 


(Y8) 


(  m  +  I ) 


sin  7TU 


sin  7 T 


m  ♦  I 


2  m  +  I 


sin 


7 T  U 
2  m+l 


sin 


77  U 
2  m  ♦  I 
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The  denominators  both  terms  rise  to  near  unity  as  u  goes 
from  0  to  m  for  most  values  of  u  they  may  be  approximated  by 
the  argument.  In  fact  for 


0  < 


tt  u 


2  m  ♦  I 


2  m  *  I 
TT  U 


sin 


TT  U 


<  1.05 


2  m  ♦  I 


2  m  +■  I 
TT  U 


Thus  we  may  approximate  the  right  side  of  (Y.8)  for  0<u<l/3  m  by 

sin  77  U  -  ■(^•{1*  sin2  (  -i-  +  — —  \  77  U  (Y9J 

7TU  77" 2  U  2  V  2  4m  *  2  / 


For  u  =  0  the  left  side  of  (Y.8)  is  at  its  maximum 


m 


(  m  +1 )  «  (  m  +  -y-) 


I  \2  (  2  m  +  l)‘ 

2 


which  is  the  limit  as  u>0  of  Y.9.  Using  (Y.8),  (Y.9)  to  approxi¬ 
mate  the  cosine  transform  of  (Y.7)  we  obtain 


.  ,  d  (  d  +  D  )  (  2  m  ♦  I ); 

Tu  =  2  '  ’u  *  *u  *  n  4 


Sin  7T  u 
IT  U 


sin 


TT  U 


7TU 


Thus  the  total  spectrum  will  deviate  from  the  average  of  its  parts  by  a 
term  whose  maximum  size  is  approximately 


m  ^ 

d  (d  +  D) 
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which  oscillates  with  a  period  of  about  2  and  which  damps  out  as 
u  increases. 

The  approximations  in  this  second  investigation  are  so 
rough  as  to  render  the  result  of  aid  only  in  determining  the  order 
of  magnitude  of  the  effect  of  this  type  of  error  on  the  behavior  of 
the  spectrum.  Since  we  will  not  smooth  by  parabolic  fit  to  segments 
in  the  future,  more  exact  work  does  not  seem  justified. 


-55- 


Section  III 

Two-dimensional  Power  Spectral  Density 
I.  Introduction 

Spectral  analysis  for  random  functions  of  one  variable 
is  common  in  engineering  work.  Problems  involving  random 
functions  of  two  variables  have  appeared  seldom  in  the  literature. 
Indeed  we  are  aware  of  but  one  article  in  which  results  of  a  two- 
dimensional  spectrum  were  presented,  "The  Directional  Spectrum 
of  a  Wind  Generated  Sea  as  Determined  from  Data  Obtained  by  the 
Stereo  Wave  Observation  Project",  Willard  J.  Pierson  cd.  (Me¬ 
teorological  Papers,  N.Y.U.  College  of  Engineering,  Vol.  2, 

No.  6,  June  1960).  The  material  in  the  present  section  is  a  further 
development  of  work  done  by  one  of  the  authors  in  connection  with 
the  above  article  and  appears  in  full  detail  for  the  first  time  here. 
The  motivations  and  interpretations  connected  with  this  work  will 
be  given  even  though  they  are  in  many  cases  repetitions  of  similar 
ideas  appearing  in  the  preceeding  sections.  Thus  the  section  will  be 
self-contained. 

To  be  able  to  speak  mathematically  about  the  height  irreg¬ 
ularities  of  a  piece  of  ground,  let  us  suppose  that  we  have  established 
a  horizontal  datum  and  a  coordinate  system  on  it.  The  coordinates 
(x,  y)  (Measurements  in  feet)  locate  a  point  on  the  datum  and  h(x,y) 
will  denote  the  height  in  feet  of  the  ground  surface  above  the  datum. 
The  pattern  of  irregularity  is  given  completely  by  h(x,y),  but,  in 
a  sense,  it  is  given  in  too  much  detail.  For  another  piece  of  ground 
having  the  same  kind  of  irregularity  the  function  h(x,y)  will  in  gen¬ 
eral  be  different.  In  order  to  extract  the  common  character  of  the 
irregularity  we  will  expand  on  the  common  mathematical  model  nf 
a  stationary  random  function. 
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The  actual  heights  above  datum  h(x,y)  include  variations 
that  are  not  of  interest--gentle  slopes,  large  hills,  etc.  We  will 
invent  a  variable  datum  called  the  "smoothed  heights"  which  con¬ 
sists  of  the  slopes  and  hills  together  with  the  general  level  of  h(x,  y) 
above  the  actual  datum.  The  differences 

d(x,y)  =  h(x,y)  ^smoothed  height  at  (x,y)j  (1.1) 

are  the  variations  of  ground  height  that  affect  the  vehicle.  We  will 
conceive  of  d(x,y)  as  being  random--that  is,  this  plot  of  ground 
gives  us  one  of  an  "ensemble”  of  possible  roughnesses  selected  at 
random  according  to  probabilities  that  characterize  the  roughness 
type.  When  we  say  that  the  general  level  has  been  taken  out  we 
mean  that  the  expected  value  or  ensemble  average,  denoted  E  (d(x,y)}, 
of  d  (x,  y)  at  any  point  is  zero,  thus,  unchanging  or  stationary  from 
point  to  point.  More  than  this  we  assume  that  all  the  characteristics 
of  the  probabilities  are  unchanging  from  point  to  point.  This  is  the 
meaning  of  "stationary  random  functions". 

Stationary  random  functions  of  one  variable  have  been  used 
in  engineering  for  more  than  20  years.  Problems  concerning  the 
output  of  a  linear  system  with  a  stationary  random  input  are  nowa¬ 
days  routine.  A  brief  summary  is  given  in  Section  II,  Part  1  of  this 
report.  In  our  work,  such  a  problem  would  be  to  determine  the  vari¬ 
able  force  on  the  driver  of  a  motorcycle  which  is  proceeding  at  a  uni¬ 
form  rate  across  our  ground  on  a  straight  line  through  (xq)  y^)  at 
an  angle  to  the  x-axis.  The  input  is  the  random  function  of  distance, 
s. 


g  (s)  =  d  (  x0+  s  cos  a  ,  yQ  a  s  sin  a  ) 


(1.2) 
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For  a  two  track  vehicle  our  input  must  be  regarded  as  two  correlated 
random  functions  (one  for  each  track);  the  formulas  are  readily- 
worked  out  as  above.  In  any  case,  the  output- -variable  force  acting 
on  the  driver  (or  the  variable  displacement  of  the  driver)- -is  also 
a  random  function.  In  engineering  practice,  the  effects  of  such  ran¬ 
dom  functions  have  been  found  to  be  related  to  the  covariance  function 
or  its  spectrum.  This  is  true  in  electrical  engineering  problems  and 
in  those  problems  of  mechanical  engineering  which  deal  with  vibration. 

In  these  problems  frequency  of  oscillation  is  of  obvious  importance. 
However,  the  spectrum  of  frequencies  has  also  been  useful  in  the 
study  of  such  "non- oscillatory"  phenomena  as  atmospheric  turbulence. 

In  this  discussion  of  random  functions  of  two  variables,  we 
will  attempt  to  use  standard  notation  and  wording  as  far  as  possible 
and  to  make  a  development  which  is  parallel  to  the  usual  one  for  functions 
of  one  variable.  This  may  be  found  in  Grenander  &  Rosenblatt,  "Sta¬ 
tistical  Analysis  of  Stationary  Time  Series"  (Wiley,  1957)  and  Blackman 
&  Tukey,  "The  Measurement  of  Power  Spectra",  (Dover,  1958). 

We  will  begin  with  the  definitions  of  covariance  and  the 
spectrum. 

R  (a  ,  b)  =  D  {  d  (x  ,y)  d  (x  +  a  ,  y  +  b)  J  (1.3) 

In  this  the  symbol  E  denotes  the  expected  value  operator  or  the  en¬ 
semble  average.  Because  of  the  stationarity  or  probabilistic  sameness 
of  one  location  to  another,  the  average  does  not  depend  on  the  location 
(x,y)  but  only  on  the  "lag"  (a,b).  The  covariance  function  is  not  ran¬ 
dom- -it  is  a  fixed  characteristic  of  the  roughness. 

a)  R  ( 0 , 0)  =  E  {(d(x,y))2}  >  0 

b)  |  R  (a ,  b)  |  <  R  (0,0) 

c)  R  (-a  ,-b)  =  R  (a.b) 


-58- 


These  properties  are  related  to  the  connection  between  covariance 
and  correlation 


COrr.  (  x  ,y) 


cov. ( x, y) 
v/var(x)  ,  var  (y) 


The  first  property,  in  other  terms,  states  that  the  variance  is  the 
same  for  any  point  x,  y.  The  second  property  is  that  of  the  corre¬ 
lation  coefficient  -  1  £  corr  (x,  y)  ^  1.  And  the  third  is  corr  (x,y)  = 
corr  (y,  x).  Since  the  correlation  is  a  measure  of  the  relationship  be¬ 
tween  the  outcomes  of  two  random  variables,  we  may  expect  R(a,b) 
to  be  large  and  positive  if  the  height  at  a  displacement  (a,  b)  from  a 
point  tends  to  be  about  the  same  as  the  heights  at  the  point  itself  R 
will  be  negative  when  large  height  at  (x,y)  is  associated  generally 
with  small  heights  at  x  +  a,  y  +  b.  The  absolute  size  of  R(a,b)  will 
indicate  the  strength  of  the  relationship.  From  continuity  of  the 
ground  height  h(x,y),  R(a,b)  should  be  near  R(0,  0)  for  small 
(a, b)-- the  correlation  between  d(x,  y)  and  immediately  adjacent 
points  should  be  close  to  unity. 

In  many  applications  the  covariance  function  has  been  use¬ 
ful  mostly  in  the  form  of  its  Fourier  transform 


F  (u.v) 


nco  -277i(xu  +  yv) 

e  R ( x , y ) d x  dy 

-oo  -oo 


(1.5) 
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Since  e  =  cos  9  +  i  sin  9,  the  quantity  2n  (xu  +yv)  must 
be  in  radians.  Since  x  and  y  are  in  feet,  2  ru  and  2 w  are  in 
radians /ft.,  or  u,  v  in  cycles /ft. 

As  with  one  variable  transform  this  can  be  inverted 

rw  x  f00  f00  -27Ti(xu+yv)  . 

R  ( x ,  y)  =  J  /  e  F(u,v)dudv  (1.6' 

-00  -CO 


The  fact  that  R  (x,  y)  is  an  even  function,  (1.4c)  implies 
that  F  (u,  v)  is  real  valued  and  in  fact 


noo 

e 

VN 


00  fCO  -27Ti(xu  +  yv) 


-00  -CO 


1/2  [  R(x,y)  +  R(-x  ,-y)  dx  dy 


.  ,  r  f00  f00  -27Ti(xu+yv) 

2  1/2  J  I  e  R (x,y)  dx  dy 

-oo  -oo 


+  y»r«t-an(-..-,.fH(Xiy)d,dy] 

-00  -00 


The  second  integral  is  obtained  by  transforming  x'  =  -x,  y'  =  -y 
then  omitting  the  primes.  The  two  may  then  be  averaged  to  give 


/•eu  /• 

(u.v)  J  J 


CO  /•  CO  -27Ti(xu  +  yv) 


1/2 


00  -CO 


R  ( x  ,y)  +  R  ( -  x  -  y)dx  dy 


,  f  r 03  C00  -  27Ti  (xu  +  yv) 

=  1/2  [  J  J  e  R(x,y)dx  dy  (1.7) 


-oo  -oo 


+ 


y0Oj.co^-27rl(-xu  -yv) 


R 


-oo  -oo 


(x.y)dxdy  J 
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which  corresponds  to  the  usual  relation  between  even  functions  and 
cosine  transforms  in  one  variable  Fourier  analysis.  From  (1.7) 
it  easily  follows  that  F(u,v)  is  an  even  function  so  that  (1.6)  may 
be  written 

,*00  /»C0 

R(x,y)  =  J  I  cos  27T  (  xu  +  yv)  F(u  ,  v)du  dv  8) 

-oo  -oo 

Some  properties  of  S(uv)  are 

a)  F  (-u,-v)  =  F  ( u ,  v) 

b)  F  (  u  ,v)  >  0  (1.9) 

/■oo  rco 

c)  J  J  F(u,v)dudv  =  R(0,0) 

-00  -00 


(Actually  (1.9a)  proved  from  the  fact  that  R(a,b)  is 

real  valued  while  (1.9b)  stems  from  a  generalization  of  (1.4a,  b) 

called  the  "positive  semi-definite"  property  of  R  to  wit,  for  any 

C.C0  C  and  points  (x., y.)  (x0,y„)...(x  ,y  )  S.  2.  C.C.R(x.  x., 
12...  n  r  '  VJY  '  2,J2'  n’Jn  l  j  xj  i-j 

■^i  ^j)$0.  The  derivation  of  (1.9b)  from  this  is  from  a  simple  ex¬ 

tension  of  Bochner's  theorem.) 
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Properties  b)  and  c)  suggest  that  the  total  variance,  R(0,  0) 
may  be  divided  up  into  components  F(u,v)  at  the  frequency  values 
(u,v).  This  is  indeed  the  case  and  one  can  derive  generalizations 
of  the  various  representations  of  a  stationary  stochastic  process  of 
one  variable.  We  will  interpret  the  spectral  density  function  F(u,v) 
as  giving  the  ”power,r  of  the  roughness  in  the  frequency  region  u 
cycles/ft.  parallel  to  the  x-axis,  and  v  cycles/ft.  parallel  to  the 
y-axis.  An  elementary  wave  at  this  frequency  will  be  a  cosine  wave 
cos  2  Ti  (ux  +vy)  which  has  the  appearance  of  corregated  iron  on  the 
xy  plane  with  the  crests  parallel  to  the  line  ux  +  vy  =  0  and  having 
crests  evenly  spaced  every  . — J  — ===—  feet.  A  representation 

V  u  *  v 

theorem  makes  us  conceive  of  d(x,  y)  as  a  sum  of  such  elementary 
waves  at  different  frequencies,  amplitudes  and  phases.  For  pur¬ 
poses  of  power  those  waves  of  the  same  frequencies  are  classified 
together.  The  power  is  analogous  to  the  sum  of  the  squared  ampli¬ 
tudes.  For  a  more  detailed  explanation,  see  Section  II,  Part  1  of  this 
report. 

A  vehicle  which  runs  over  the  rough  ground  is  affected  only 
by  the  roughness  under  its  wheels.  If  its  track  is  straight  across  a 
stationary  field,  the  roughness  is  stationary.  The  simplest  case 
deals  only  with  one  track  but  more  realistically  we  will  have  two 
parallel  tracks.  Let  us  derive  first  the  formulas  for  the  single  track. 
Equation  (1.2)  gives  the  random  function  involved.  The  covariance 
function  for  this  is 


C(r)  =  E  {  g(s)g(  s+  r)} 

=  e  jd  (x0+  s  cos  a,y0+  s  sin  a)  d  (  x0  +  scos  a  +  r  cos  a  , 

y0+ s  sin  a  +  r  sin  a  )|  (i-i°) 


=  R  ( r  ccs  a  ,  r  sin  a) 
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The  spectral  function  for  C  (  t)  is  the  function  D  (  X)  that  satisfies 

C(r)=/V,nXrD(\)d\ 

-00 


But  we  have  from  (1.6) 


R  ( r  cos  a  ,  r  sinajr/*00  f°°e 
-oo  -oo 


27Ti(ur  cos  o  +  vsin  a) 


F  (u,v)  du  dv 


and  if  we  take 

u  cos  CL  +  v  sin  a  =  X 

-  u  sin  a  +  v  cos  a  =  X 


(1.11) 


whence 

X  cos  a  -  u  sin  a  =  u 
X  sin  a  +  u  sin  a  =  v 

and 


J 


d(X,u) 

cHu,v) 
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the  integral  may  be  transformed  to 


f°°e27ri^r  I"  J  f(Xcos  a  -  fJ.  sin  a  ,  X  sin  a  +  fj.  cos  a)du  ldX 

-CO  -00  J 


from  which  we  have 


D(X)  =  J  F(X  cos  a  -  fJL  sin  a  ,  X  sin  a  +  /x  cos  a )  du 


(1.12) 


relating  the  linear  spectra  in  various  directions  to  the  two  dimen¬ 
sional  spectrum. 

For  parallel  tracks  separated  by  a  ■//P\x  y  ) 

distance  b  the  spectrum  of  each  track  is  y//'' 

given  by  (1.12).  The  cross  covariance  of 
the  two  profiles  is 

H(r)  =  E  {d(x0+  s  cosa,  yp+sin  a  )  d  ( x,  +  s  cos  a  +  r  cos  a ,  y,  +  sina+  r  sin  a)} 
=  R  (  x,-  x0  +  r  cos  a  ,  yt  -  yQ  +  r  sin  a) 


If  (x^)  is  opposite  (xoyo)  011  a  line  perpendicular  to  the  track 
x,-x  =  -b  sin  a  y,-y  =  b  cos  a.  Since  H(t)  is  not  an  even  func- 
tion,  its  Fourier  transform  has  both  a  real  and  an  imaginary  part. 
These  are,  respectively,  its  cosine  and  its  sine  transform.  They 
give  the  frequency- -phase  difference  relationship  between  the  two 
tracks.  As  above,  this  complex  spectrum  is  the  function  K  {X) 


such  that 


r oo  27nXr  .  , 

H (r) J  e  K (X) dX 
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and  (1.12)  and  (1.6)  gives  us 


l_l  ^  j  J*C°  yO0^277"i  [("bsinft  +  rcosa)  u  +  (b  cos  a  +  Tsina)  v]F(u,v)du  dv 


-00  -00 


Making  use  again  of  the  transformation  (1.11) 


/•co  /*oo  27Ti  ( bu  +  rx) 

H (r )=/  J  e  F(Xcosa-^sina,X  sin  a  +/xcosa)du  dv 

-oo  -oo 


whence 

/•co  2rrib/i 

K(X)=  J  e  F  (cos  a sin  a  ,  X  sin  a  +  fJ-  cos  a)  du 

-co 


Formula  (1.12)  is  a  special  case  of  this  where  b  =  0.  For  each  X 
the  integral  is  a  line  integral  of  the  two  dimensional  spectrum  along 
lines  perpendicular  to  the  direction  (a)  of  the  track  and  at  a  distance 
from  the  origin. 

This  work  gives  two  possible  derivations  that  may  be  made 
from  the  two  dimensional  spectrum.  In  both,  the  track  of  the  vehicle 
is  straight.  When  the  track  curves,  the  resulting  roughness  is  not 
stationary  unless  the  roughness  of  the  field  is  the  same  in  all  direc¬ 
tions.  These  spectra  (1.14)  for  various  a  will  give  the  conditions 
met  by  the  vehicle  in  various  directions  and  may  be  of  aid  in  study¬ 
ing  the  effects  of  slowly  varying  non-  stationarity. 

In  this  introduction  we  have  dealt  with  the  theoretical  quanti¬ 
ties  associated  with  the  stochastic  process.  We  have  described  them 
and  indicated  how  they  can  be  used.  Our  present  problems  are  more 
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concerned  with  estimating  these  quantities,  particularly  the 
spectral  function.  Estimation  of  a  two-dimensional  spectrum  is  a 
simple  extension  of  one  dimensional  spectral  estimation.  Since 
the  formulas  are  not  available  in  the  literature  we  will  present  a 
fairly  complete  account  of  them  here. 


2.  Relation  of  the  Theory  to  the  Practical  Problem. 

As  in  all  statistical  work  we  will  use  data  from  outcomes 
of  the  random  scheme  to  make  estimates  of  the  desired  character¬ 


istics  of  the  probabilities.  In  our  case  the  desired  characteristic 
is  the  spectral  density  function  of  the  covariance.  Statistical  estima¬ 
tion  makes  use  of  the  fact  that  when  repeated  outcomes  are  averaged, 
the  random  part  of  the  average  tends  to  be  small.  This  repetition 
may  be  accomplished  with  d(x,y)  by  taking  the  values  at  different 
places;  the  stationary  property  is,  in  effect,  the  property  that  the 
randomness  repeats  itself  in  different  places.  Actually  a  some¬ 
what  more  stringent  property,  the  ergodic  property,  is  necessary 
so  that  averages  over  (x,  y)  tend  to  ensemble  averages.  We  will 
assume  as  do  most  practitioners  that  this  property  holds. 

It  is  most  convenient  for  computation  and  for  measurement 
to  take  data  at  intersection  points  of  a  rectangular  lattice.  Though 
the  data  we  have  taken  has  the  separations  equal  in  both  directions, 
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Supposing  the  smoothing  of  the  profile  to  have  been  done,  a  topic 

we  will  discuss  in  Section  5,  our  data  will  consist  of  n  by  n 

x  J  y 

values 

d  (k  Ax  ,  j  Ay) 

Only  a  part  of  the  randomness  is  involved  in  this  data.  The 
covariances  which  underlie  it  are  of  the  form 

E  |d(k  Ax  ,  j  Ay)  d((k  +  a)  Ax ,( j  +  b)  Ay)}  =  R(a  Ax ,  bAy)  (2.1) 

for 

a  =  0,  ±  I ....  ±  nx  b  =  0,  ±  I , . . .  ±  ny 

No  other  values  of  the  covariance  function  are  used  to  "pro¬ 
duce"  the  data  and  we  cannot,  therefore,  hope  to  make  direct  esti¬ 
mates  of  any  but  these. 

Supposing  that  we  knew  these  values  of  R(x,y),  what  can  we 
find  out  about  the  spectrum  ?  This  is  a  preliminary  question  and  we 
must  later  refine  the  answer  because  we  will  not  have  the  values  of 
R  (x,y)  but  only  estimates  of  them.  This,  of  course,  will  affect  the 
way  we  answer  the  preliminary  question. 

The  spectrum  is  related  to  these  values  of  R(x,y)  by  equation 

/CO  rCD 

J  cos  27T  (a  Axu  +  b  Ayv)  F(u,v)du  dv  (1>8) 
-co  -co 


(1.8) 


-67- 


In  analogy  to  equation  (1.7)  we  will  later  multiply  the  R(a  A  ,  bA  ) 

x  y 

by  cosines  and  add.  But  we  want  to  investigate  the  matter  generally 
so  instead  of  using  cosines  to  obtain 

Z  Z  cos  27T(q  Axu0  +bAyV0)  R  (a  Ax ,  bAy)  Ax  Ay  (2  0) 

let  us  use  general  coefficients  to  obtain 

F*(u0,v0)  =  £  £  w  (a  ,  b  u0, v0)  R(aAxbAy) 
a  b 

=  Z  I>(a,b  ;u0,v0)  J  J cos  2  7T  (a  Axu  +  bAxv)F(u,v)du  dv 
a  b  J  J 

(2.3) 

■//{?£«<  a,b  ;  u0,v0)  cos27r(a  Axu  +  bAxv)  F(u  ,v)du  dv 


From  this  we  see  that  F*  (uo,vq)  is  a  weighted  average  of  spectral 
values, 

F  *(u0,v0)  =  J"Jw(u0,v0  ;  u,v)  F(u,v)  du  dv  (2.4) 

Our  purpose  is  to  obtain  F(uq,v  ),  therefore 

W(u0,v0  ;u,v)  =  £  £w(a,b;  u0  ,v0)  cos  27T  (a  Axu  +  b  Axv  ) 
a  n 


(2.5) 
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is  ideally  a  Dirac  or  delta  function  of  (u,v)  centered  at  (uq,  v  ). 

It  is,  of  course,  impossible  to  attain  this  ideal  since  W(‘)  is  a 
finite  Fourier  series.  An  obvious  W(‘)  to  try  is  the  finite  part  of 
the  Fourier  series  for  the  Dirac  function 


n 


x 


I 


a  =  -  nx 


Ax  Ay  cos  27r(aAxu0  +  b  Ayv0)  cos  27r  (a  Axu  +  bAyv) 


(2.6) 


which  corresponds  to  (2.2). 

We  will  not  go  into  detail  here  about  the  form  of  W(uqvo;u,v). 

Let  us  note,  however,  from  (2.5)  that  it  is  even  and  periodic  with  a 

period  of  1/A^inthe  u  direction  and  1/A^  in  the  v  direction  so 

that  if  it  has  a  sharp  peak  at  (u  v  )  it  has  identical  peaks  at  (+u  +r/A  , 

+  vq  +  s/A^)  for  integers  r,s.  Referring  to  formula  (2.4)  we  sec 

that  the  sharp  peak  of  W(*)  at  (uo,vq)  picks  up  the  values  of  F(u,v) 

near  (UQ.V0)  while  the  other  peaks  get  the  values  near  their  positions. 

F*(u0,v  )  is  therefore  made  up  of  an  average  of  the  values  of  F(u,  v) 

near  (u  ,  v  )  with  the  addition  of  similar  averages  for  values  near  the 
o  o 

other  peaks  of  W(').  These  other  averages  are  called  "aliases". 

Their  presence  is,  of  course,  unwanted.  It  is  due  to  the  fact  that  the 
data  is  taken  at  equally  spaced  points.  The  smaller  the  spacing 
(A^,  A^  the  farther  apart  are  the  peaks  of  W(’).  The  usual  spectrum 
"dies  out"  for  large  (u,v)  and  if  A^,A^  are  small  enough,  the  aliased 
values  of  F(u,v)  will  be  at  large  enough  frequencies  that  they  do  not 
affect  the  estimates.  The  periodicity  of  W(’)  means  that  the  function 
F*(u,v)  of  (2.3)  is  itself  periodic.  Actually  the  values  outside  the  rec¬ 


tangle 


(- 


■zrr-  <  u  <  — -7- 

2AX  -  -  2 A, 


<  v  < 


2  Ay  ~  _  2  Ay 


(2.7) 
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are  repetitions  of  those  inside  this  rectangle.  Moreover,  the  spec¬ 
tral  function  is  even  (1.9a).  Consequently  all  the  values  of  the  es¬ 
timated  spectrum  may  be  found  from  those  within  the  rectangle 

(-1/2A  ^  u  $  1/2A  ;  0  v^  1/2A  )  and  these  need  be  the  only 
xx  y 

values  reported.  (For  a  one-variable  spectral  estimate  the  only 
values  ordinarily  presented  are  for  positive  frequencies  out  to  the 
Nyquist  frequency  of  1/2A. ) 

The  phenomenon  of  aliasing  is  a  limitation  on  accuracy  re¬ 
sulting  from  sampling  at  intervals.  It  arises  fromthe  form  of 
Eqn.  (2.5)  and  not  from  the  specific  coefficients  w(a,b;  u  ,v  ). 

Another  limitation  of  accuracy  arises  from'  our  having  only  a  finite  number 
of  data  points.  Equation  (2.5)  is  a  Fourier  expansion  with  all  co- 

Since  the  coefficients  may  be  cal- 


efficients  zero  for  a$n  b$n 

^  y 

culated  from  W(‘ )  by 


w  ( a  ,  b,u0  v0)  =  Ax  Ay 


/:- 1 


2  A, 


(2.8) 


2  A, 


2  A. 


w(u0,v0;u  v)  cos  27T(a  Axu  +  bAyv)  du  dv 


it  is  clear  that  the  ideal  form  for  W(‘  )  ,  (Dirac  "delta"  function), 

cannot  be  realized,  otherwise  all  w  (  )  would  be  non-zero.  In 

fact  if  W  (  )  is  to  be  large  near  both  positive  and  negative  parts 

of  cos  2  7r  (a  A  u  +bA  v)  for  a>n  ,  b^n  for  the  integrals  to  be 
xx  x  y 

zero.  This  indicates  that  the  nori-zero  area  of  w(‘  )  should  be  at  least 
a  half  period  wide  and  a  half  period  high  or 


_J _ 

2  nx  Ax 


I 


by 


2ny  Ay 


(2.9) 
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As  a  cheek  on  this  we  note  that  the  top  line  of  (2.3)  connects  a 
finite  number  of  R  values,  by  a  linear  relation,  to  the  values  of 
F*.  Under  any  but  abnormal  conditions  an  equal  number  of  values 
of  F#  determine  all  the  R's  and  hence  also  the  other  values  of 
F*.  If  we  space  these  values  evenly,  in  the  relevant  area  of  fre¬ 
quencies,  (-  1/2A  £  u  $  1/2A  ,  0  $  v  <  V  2A  ),  the  spacing  will 

x  x  y 

be  l/2n  A  by  l/2n  A  as  before, 
x  x  J  '■  y  y 


3.  Statistical  Problem 

In  Section  2  we  acted  as  though  the  values  of  R(aAv,  bA 

were  known  for  a  =  0.  +  1.  .  .  +  n  b  =  0,  +  1 

—  x 


y> 


+  n^.  Actually, 

we  will  estimate  these  values  from  the  height  data  dfkA^,  jA^).  The 
usual  estimate  for  an  expeeted  value  (2.1)  is  an  average 


(nx- la  I )  (ny- |  b|)  Z  Z  d  ( k  Ax ,  jAy)  d(k  +  a)  Ax  ,(j+b)Ay  =  ra,b  ^3- ^ 


where  k  and  j  range  over  all  the  pairs  on  hand.  The  aetual  limits 
of  summation  depend  on  the  signs  of  a  and  b.  Under  the  ergodie  hy¬ 
pothesis  this  converges,  as  nx  and  n^  grow  larger,  to  RfaA^, 

bA  )  For  values  of  a,  b  close  to  n  ,  n  ,  the  sum  has  only  a 

y  x  y  J 

few  terms  and  the  estimate  is  not  good.  Therefore  we  must  be  eon- 

tent  with  estimates  r  .  for  -m  ^  a  <  m  ,  -m  ^  b  <C  m  where 

a,  b  x  v  x  y  y 

m  and  m  are  considerably  smaller  than  n  and  n  .  The  expected 
x  y  J  x  y 

value  of  r  ,  is  R(aA  ,  bA  ).  The  different  r  .  are  averages 
3. ,  D  x  y  a,D 

of  different  numbers  of  terms,  those  of  larger  a  and  b  have  fewer 
terms  and  thus  less  precision.  In  general,  also,  the  different  r's 
are  correlated. 
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In  accordance  with  our  findings  in  Section  2,  we  will  try 
to  estimate  values  of  the  spectral  function  which  are  evenly  spaced 
throughout  the  relevant  range.  Since  we  will  have  covariance  values 
only  up  to  m  ,  m  ,  the  spacing  of  spectral  values  is  correspondingly 
larger.  For  mathematical  convenience  we  will  use  a  spacing  of 
l/(2m+  1)A  rather  than  l/2mA,  where  m  is  in  the  x-direction 
and  my  in  the  y- direction.  Thus  we  will  try  to  estimate  (see 
Eqn.  (2.3) 


a 


A. 


(2  mx+  I )  (2  my  +  l ) 


=  F* 
a>/3 


(3.2) 


Q  =  0,  ±  I,  ±  2 . ± 


+  rrtw 


j3=  0,±l  ,  ±  2 . ±  my 


Using  the  coefficients  w(a,b;  uo>  v  )  of  Eqn.  (2.2)  we  can  calculate 
the  estimates  ... 


a./3 


m*  m 

'  I  I 

a=-mx  b=-mv 


w 


(a.bi - 9L 

'  (2m.+  l 


_  ,  P 

l)Ax  (  2  my+  |)  Ay  )  rQ,b  (3.  3) 


whose  expected  values  are  F*  0. 

a,  P 

To  approximate  the  distribution  of  the  random  variable  f^  ^ 
we  will  assume  the  Gamma  distribution  which  has  been  found  to 
fit  well  in  this  kind  of  work.  The  parameters  A,K  of  the  distribu¬ 
tion  will  be  estimated  by  the  method  of  moments. 
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The  following  is  the  mathematical  expression  of  our  first 


assumption: 


„  f  .  .  )  AK  /•  A  K-l  -  At 

P  \  f„  a  <  A }  =  -=r —  /  t  e 

i  j  f(K)  J0 


As  in  Section  II,  Part  3  above,  we  have  corresponding  to 


(3.2) 


(3.4) 


Confidence  limits  may  be  found  as  in  Section  II,  Part  3  and  inter¬ 
preted  by  means  of  the  table  there. 

The  derivation  of  the  formula  for  Var  {f^  ^  which  also  gives 
cov  (f^ p  p}  =0  is  somewhat  tedious  and  will  be  presented  in  an 
appendix  to  this  section.  Several  simplifying  assumptions  are  necessary 
and  since  these  are  relevant  in  practice  we  will  discuss  them  below. 

The  result  is  a  formula  for  the  degrees  of  freedom. 


deg.  fdm.  =  2  K  «  2  (  j 


(3.  5) 
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where  g  =  m  /n  =  m  In  is  the  ratio  of  the  number  of  R  values 

xx  y  y 

in  each  direction  the  number  of  profile  values  in  each  direction. 

The  assumptions  made  in  the  approximate  calculation  of 

var  {f  0  }  and  cov  {  f  0,f  .  0.  }  are  as  follows 
a,  p  J  a,  p  a' ,  p'  1 

A. 2  The  profile  height  deviations  d(x,y)  of  (1.1)  are  approximately 
Gaussian  with  means  0. 


A. 3  The  values  of  m  and  m  are  small  in  comparison  to  n 

x  y  x 

and  n  so  that  the  numbers  of  terms  in  the  sums  for  different  co- 

y 

variance  estimates  (3.1)  are  approximately  the  same. 


A. 4  F(u,v)  is  essentially  zero  for  (u,  v)  outside  the  range  (2.8) 
so  that  aliasing  does  not  occur.  It  is  relatively  constant  within 
any  of  the  small  intervals  (2.10). 


A. 5  The  non-zero  parts  of  the  functions  W(‘)  of  Eqn.  (2.4)  are  no 

wider  than  (-^  ,  by  (-  ),  so  that  their 

xx  xTx  y  y  y  y 

"overlap"  is  negligible. 

We  have  little  control  over  assumptions  A.l  and  A. 2.  As¬ 
sumption  A. 5  is  somewhat  in  conflict  with  the  discussion  of  Eqn.  (2.  9) 
which  says  that  the  non-zero  part  of  W(‘  )  must  be  at  least  that  wide. 
To  come  near  a  compromise  between  A. 5  and  this  discussion  wc 
must  use  a  relatively  long  Fourier  expansion  (2.5)  and  choose  the 

coefficients  with  care.  This  means  a  fairly  large  m  ,  m 

x  y 

also  desirable  for  good  resolution  or  closeness  of  spectral  esti¬ 
mates.  But  A. 3  requires  an  even  larger  (n  ,  n  ).  The  ratio  g,  of 

X  j 

m  to  n  should  also  be  small  to  obtain  a  large  degree  of  freedom 


(3.. 


Assumption  A. 4  relates  to  the  size  of  the  spacings  A^.,  Ar 


in  the  profile  data. 
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The  desire  for  good  resolution,  statistical  precision,  and 
a  wide  range  for  the  values  of  F(u,v)  require  a  large  number  of 
profile  values  that  are  closely  spaced.  The  expense  of  measuring 
the  profile  and  handling  the  data  must  be  balanced  with  these  desires. 
Even  if  large  amounts  of  data  could  be  obtained  and  processed  cheaply, 
the  lack  of  stationarity  and  the  inherent  errors  of  measurement, 
which  have  not  been  taken  into  account  here,  would  limit  the  accuracy 
and  precision  of  the  results  in  other  ways  so  we  will  always  be 
faced  with  the  above  dilemma.  The  principles  outlined  in  this  sec¬ 
tion  will  help  us  to  deal  rationally  with  it. 
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Appendix  to  Section  3 

Derivation  of  degrees  of  freedom  formula. 

Let  us  consider  two  quantities  of  the  form  CL  3) 
m.  m„ 


f.  ■ 


i  2  ^ 

a  =  -mx  b  =  -  my 
mx  my 

V  ^  T  hQb  'ob 

a  =- mx  b=-my 


cov 


{f|f2}  =  I  I  I  I  SQb  htfb.  cov  (rQb  ,  rQ.b,)  (xi) 


In  this  sum  the  terms  cov  (ra^>  ra,.  ,)  are  readily  worked 
out  using  (3.1) 


cov  ( ra  b  ra'b' )  = 


(nx-  I  a  I)  ( ny- 1  bl)  (nx-  |a'| )  (n  -  ib'l) 


Ij  Ij  Z,i  Ij.  CCi.j 


(X2) 


where  the  limits  of  the  sums  depend  on  a,  b,  a1,  b'  and 

C(i,j  ;i\j')  =  cov  |d( i Ax  ,  j  Ay)d  [(i  +  a)  Ax,(i  +  b)  Ay  j  , 


d  ( 


i'Ax  j' Ay)  •  d  [(i'  +  a1)  Ax  ,  ( j'  +  b’)  Ay  H 


(X3) 
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Using  the  Gaussian  assumption  A. 2  we  can  evaluate  this 
covariance  from  the  formula 


cov  (MN.ST) 


=  E  {  [  M  N  -  E  ( M  N)  ]  [ST-E(ST)]} 


,  mnst} 

-  E  J 

[mn] 

•  e(st} 

=  E  • 

{ms}  e 

M 

>  +  E 

{mT}e{ 

In  this  the  Gaussian  assumption  is  used  only  to  evaluate  the  fourth 
moment. 

Using  this,  (X.3)  becomes 

C  (i  j  ;i‘ j')=  R  [  (i'-i)Ax  ,(j’-j)  {  (i'-i  +  a’-a)  AK  .(j'-j  +  b1- b)  Ayj 

(X4) 

+  R  [  ( i’+a'-i)  Ax,(j'+b’-j)Ay  ]  R  [  ( i'-i -a)  Ax ,( j‘- j -b Ay)  ] 


We  need,  as  our  result,  (X.l)  in  terms  of  F,  so  we  will 
apply  (1.6)  to  express  the  product  of  R's  as  an  integral 


r  r  r  r  27ri(ou«-  bv  ♦a'u'  +  b'v1) 
R  (ab)  R(a'b‘)=/  j  J  J  e 


F  (uv)  F  (u1  v’ )  du  dv  du’  dv' 


Using  (X.5)  in  (X.4)  we  obtain 


(X5) 


C(ij;i'j')=  J  J  J*J*EF(uv)F(u'v')dudvdu’dv’ 


(X6) 
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where 

E  =  exp  27ri  [  ( i'-  i )  Ax  u  +  ( j'-  j  )Ay  v  +  (i,-i+a,-a)Axu'+(j'-j  +  b'-b)  Ay  v’  ] 
=  exp  2  7ri  £  ( i'+a'-i )  Axu  +  (j'+b'-j)  Ay  v+(i'-i-a  )  Axu'+(j'-  j-a  Ay)v'j 

=  exp27ri  £  ( i‘  -  i )  Ax(u  +  u' )  +  ( j'- j )  Ay  (  v  +  v' ) -a  Ax  u'  -  b  Ay  v'  j  x 
x|exp27ri  ^  a1  Ax u'  +  b' Ay v'  j  +  exp  2 7ri  ^a'Axu+b'Ayvj| 

=  2exp27Ti  (i‘-i+  )Ax(u'+u)+(j'-j+  -|-)(v,+  v)-a  Axu'-b  Ay  v '  1 

cos  27T  (a'Ax  U  g  u~  +  b'  Ay  v  2V  ) 

Because  F(u,v)  is  an  even  function,  the  integral  (X.6) 
with  u  interchanged  with  v  and  u'  with  v'  also  equals  C(i,  j;i\  j’). 

It  has  in  place  of  E,  a  function  identical  with  E  except  for  the 
last  two  terms  of  the  exponential.  The  average  of  the  two  integrals 
also  equals  C(i,j;i',j')  and  is  like  (X.6)  except  that  in  place  of  E 
it  has 


exp  27Ti  ■[(i'-i+  Ax(u'+u)  +  (j,-j  +  -|-)Ay(v'+v)}cos  27r(a’Ax  U  2U  +  b'Ay 
Jexp27Ti(-a  Axu'-bAyV' )  +  exp  2  7n(-aAxu-bAyv)  J 

(X  7) 

=  2exp27Ti  [(i'-i+  Ax(u'  +  u)+(j'-j  +  b  gb-  )  Ay  (v'  +  v)j 

cos 2rr ( o  Ax  -u-|-u-  +  b  Av  )cos27t (a’A*  f  b'Av  ■v-2' ~ ) 


iy  2 


x  2 


iy  2 
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We  will  imagine  this  substituted  into  (X.6)  and  thence 
into  (X.2)  where  we  must  sum  the  integrals.  The  sums  are  finite 
and  can  be  taken  into  the  integral  where  they  are  found  to  affect 
only  the  exponential  factor  of  (X.7).  This  sum  falls  into  the  product 
of  four  sums  of  the  form 


27T  i  A  x(  u'  ♦  u)  i1 

I  e 

i' 


(X  8) 


and  a  factor  ( 

2171  [■^  AX(U'*U)+  -^-Aytv'  +  V)] 

e 


(X  9) 


The  sum  (X.8)  is  a  geometrical  series  and  its  limits  depend  upon 
the  value  of  a.  It  can  be  put  into  the  closed  form 

27Ti  Ax(u  +  u')i  27Ti  Ax  (nx-o'  +  l) 

e  =  e 


sin 27T Ax  (nx-  I  a  |) 

sin27rAxiLyi:- 


The  other  sums  may  be  worked  out  similarly.  Their  product  con¬ 
sists  of  four  factors  that  are  sine  ratios  and  an  exponential  factor 
which  will  be  found  to  be  the  reciprocal  of  (X.9).  Thus,  the  formula 
(X.2)  for  the  covariance  of  two  7  values  becomes 


_ 2 _ 

(  nx-  |  a  |  )  ( nx  -  |  a'|  )  (  ny  -  1  b  | )(  ny  -  j  b’l ) 


//// 


sin  7T Ax(u  »  u*)  (nx-  1  a1 1)  sin7T Ax(m-u')  (nx-  1  0  j)  (X  10) 


sin7T  Ax  (u  ♦  u‘) 


sin7rAx(u  +  u’) 
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sin  7T  Ay  ( v  ♦  v1)  (uy  -  1  b1 1)  sin7rAy  (v»v|)(Uy-  |  b|) 

sip  7T  Ay (v ♦  V1 )  Sin  7rAy(v+v')  (X  10  con.  ) 

cos27r(aAx^  +  bAy 

cos  27T(a'  Ax  ■  u-g-—  +  b  Ay  -■  y-  F  (u,v)  F  (u1  v1)  du  dv  du1  dv1 


Before  substituting  this  long  expression  into  (X.l)  where  it  must 
be  multiplied  by  coefficients  and  summed,  let  us  examine  the  in¬ 
tegrand.  The  sine  ratios  may  be  combined  with  the  first  factor 
to  give  four  functions  of  the  form 

sin  (a£E) 

N  sin  A  £ 

Each  of  these  has  a  spike  at  £  =  0  of  height  1.  and  also  spikes 
at  §  ~  -  for  the  integer  K.  The  factor  N  affects  only  the  width 

of  the  spikes.  If  we  adopt  assumption  A. 3  the  dependence  of  the 
sine  ratio  terms  on  a,b,a',b'  largely  disappears.  The  product 
of  four  factors  is  approximately  a  product  of  two  squares,  and  the 
indices  of  summation  in  (X.l)  affect  only  the  cosine  terms  of  (X.10). 
Let  us  take  the  sum  and  the  factors  g^  h^^,  inside  the  integral 
whence  we  have  the  sum 


z  z 

o  b 


gab  cos  Z1T  (a  A, 


u'-  u 


+  b  A, 


V'- V 


)Z  Z  h  ,b,  c03  27r(clA^u-^y-+b'Ay1^-) 

a'  b'  d  y  2 


=  G  ( 


u'-u 

2 


2 


2 


) 
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where  G  and  H  are  two  of  the  functions  W  defined  in  (2.4), 
i.e.,  W  for  two  different  pairs  (uovq). 

We  may  now  express  (X.l)  as  an  integral  of  the  spectral 
function.  Let  us  denote  the  "common  value"  of 


nx  -  I  o  |  =  nx  -  |  o'  |  =  nx 
ny  -  1  b|  =  ny-  |  b'|  =  Ny 


cov(f,f2)=  2  ////[ 


sinTrA^u+u'JN, 
Nx  sin7T  A.(u+u 


nsin7rAy  (v  +  v')nv  y 
Ny sin7T Av(v  +  v')  J 


G  (  »  ,V-  2~— )  H  (  u  2U  *  -~-2  )  F  (u  v)  F(u'  v1)  du  dv  du'  dv1 


Let  us  change  variables  so 


“(u'-u) 
1 


r 


•j(v'*v) 

_L(v'-v) 


d  ( w  w1,  z  z1 ) 
d  (u  u1 .  v  v1 ) 


=  |J|= 


=  w 


=  w' 


=  z 


-  z' 


T  i  00 
00 

00 

00  ±-4r 


w  +w' 
w  -w' 
z  +  z' 
z  -z1 


4 


_l_ 

4 


=  v' 


=  V 


,x  <  .  «  CCCCt  sin£7TAxNxw  sz ,  sinE7rAvNvz 

C0»«,.f,).2 JJJJ  (-^;s27rA<w  )  ( 1T-Sin JN;Z  )  St*M'IHI»',!') 


X  F(w- w1 ,  z-z')  F  (  w  +  w' ,  z+z1)  4  dw  dw'  dz  dz' 
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Carrying  out  the  integration  first  with  respect  to  w  and 
z  we  note  that  the  squared  sine  ratios  have  spikes  for  w  (or  z) 
equal  to  0,  _+_  (or  0,  +  ■ )  etc.  and  that  except  for  these 


x  yl  1 

spikes  whose  width  is  about  2*  a  (  or  2  a  )  the  values 

XX  y  y 

are  quite  small.  Using  assumption  A. 4,  F  is  relatively  constant 
inside  the  peak  at  0  and  has  very  small  values  in  the  region  of  the 
other  peaks.  The  first  two  integrations 


G(w'z')  H  (w'z1) 


[f(w'z')] 


sm  27T  Ax  N,w  \2  , 

— — : - —  )  dw 

Nxsin27rAxw  ' 


The  two  integrals  of  the  sine  ratio  may  be  evaluated  from 
Fejer's  integral  (Titschmarsh  "Theory  of  Functions",  p.  413).  The 
result  is 


2 


f  f  G  (w1  z 


)  [ F ( w 1  z 1 )  1  dw1  dz1 


dw'  dz1 


COV  (  f,  ,  f  2  )  = 


Ax  Ay  Nx  Ny 


(X  12) 
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From  assumption  A. 5  when  G  and  H  are  two  different 
w  functions  their  non  zero  parts  do  not  overlap  and  so  the  co- 
variance  of  two  spectral  estimates  is  zero. 

The  variance  of  f  is  the  covariance  of  f  with  itself  and 
formula  (X.ll) 


var  ( f )  = 


2 

Ax  Ay  Nx  Ny 


//[  G(w'  z1)  F(w'z 


>] 


2 

dv1  dz1 


Let  us  recall  that  our  purpose  is  to  compute  the  degrees 
of  freedom  using  formula  (3.5).  Our  result  (X.12)  will  be  in  the 
denominator  and  the  numerator  will  have 


[*{*}]■[//  G  {w\  z1)  F(w\  z')  dw'  dz'  ] 


(X  14) 


This  follows  from  taking  the  expected  value  on  both  sides  of  (3.3) 

then  using  (2.2).  In  both  (X.12)  and  (X.13)  we  will  take  G(wz') 

to  be  an  abbreviation  of  w(uo>vq;  w'.z')  of  (2.4) 

To  attain  approximations  of  the  two  integrals  we  will  make 

strong  use  of  assumption  A. 5  and  also  assume  an  idealized  rectangular 

form  for  G.  ’  This  form  proceeds  from  our  discussion  at  the  end  of 

Section  2;  it  is  that  G  is  0  in  the  relevant  region  of  frequencies 

(w',z')  except  in  two  rectangles  of  dimensions  — ^ -  X  — ^ - 

x  x  y  y 

centered  at  (u  v  )  and  (-u  -v  ).  The  total  volume  of  these 
o  o  oo  2 

rectangles  should  be  unity  so  the  height  of  G  will  be  -5  (m  m  A  A  ). 

“  x  y  x  y 

With  this  (X.12)  and  (X.13)  work  out  to  be  approximately 


var 


(’}  *  ~/&:2n;n.  [^“ov,,)]2  a,  a,) 


*X  y  *’X  "y 
1  2 


(X  15) 


[  E  {  f  }  ]  =  [  F^o  Vo)] 


(X  13) 
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The  degrees  of  freedom  formula  (3.5) 

r  iZ 

2  I  F(u.  v.l 

d.  f.  =  2K 


[f  <“.*.)] 


2  m »  m  v 

NxNy 


Nx  Nv 

Mx  My 


Now  if  we  recall  the  definition  (X.ll)  of  Nx 
an  approximation 


and 


we  have  as 


d.f. 


n x  mx  ny  ~  ^y 
mx  my 


which,  when  we  use  the  fixed  fraction  g  for  both  m  and  m 

B  x  y 

we  get  the  formula  (3.6).  This  approximation  is  probably  as  good 
as  the  assumptions  justify  but  more  exact  results  may  be  found 
by  approximating  more  closely  (X.13)  and  (X.14)  using  the  actual 
W  function  instead  of  the  rectangular  idealization. 
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4.  Filtering 

The  assumptions  at  the  end  of  Section  3  are  used  in  both 
the  statistical  work  and,  as  we  pointed  out  in  the  discussion  of  them, 
in  relating  the  theoretical  work  to  the  practical.  We  will  discuss 
here  assumptions  A. 4  and  A. 5  beginning  with  the  latter. 

The  shape  of  the  function  WO)  should  be  two  symmetri¬ 
cally  placed  peaks  on  a  field  of  zeros.  Though  considerations  in 
Section  2  indicat  that  the  peaks  should  be  as  narrow  and  high  as 
possible  around  the  central  point,  the  statistical  considerations 
expounded  in  equations  (X.14)  and  (X.15)  indicate  that  the  variance 
of  f  increases  with  the  narrowness  of  the  peak.  Thus  we  are  led 
to  seek  a  box  shaped  W(')  function.  The  choice  of  this  function 
is  the  subject  of  a  good  deal  of  literature  .  Here  we  will  follow 
some  of  the  early  work- -that  given  in  Blackman  and  Tukey, 
"Measurement  of  Power  Spectra".  This  will  take  the  direction  of 
using  as  standard  the  W(*)  given  by  (2.6)  and  modifying  it  by 
filtering.  Let  us  note  first  several  formulas  concerning  finite 
Fourier  expansions. 

Suppose  we  have  an  expansion  C(u,  v)  with  coefficients 

c  , 
a,b 


C  (u  ,  v  ) 


my 

I 

b  =  —  my 


27T  i  (a  Ax  u*b  Ay  v ) 


'ob 
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and  we  want  to  select  evenly  spaced  values  of  (u,v)  to  "represent" 
W(u,  v).  Let  us  take 

u  =  0  .  ±  8X  .  -  2BX  ,  . . •  ,  ±  mx  8X 

v  =  0  i  i  Sy  »  —  2§y  > .  .  •  * i 


m y  Sy 


rn^  my 

Cj,k  =  C(  jSx,  kSy)  =  A  I  cQ 

o=-mx  b=-mv  ' 


27H(  AXSX  Oj+AySyby) 


Multiplying  both  sides  by 

27Ti  (AXSX  +  AySy  bk) 
e 


and  summing  on 


(i.k) 


gives 


27Ti  (  AXSX  aj  +  Ay  Sy  b  k  ) 

Cj.k  e 


=  X  l 


'a,b 


27Ti(AxSx  j  (a -a  )+  AySy  k(b-b)) 
e 


-86- 


The  second  sums  are  easy  to  work  out.  The  exponential  is  readily 
separated  into  a  product  and  gives  a  product  of  two  geometric 
series,  each  of  the  form 

m  ASi  277-iAxSxmx(o-a)  27Ti  Ax8x(mx  +  l)(a-a ) 

mx  27Ti  AxOx  Ha -a )  e  -e 

.  Z  e  ‘  ’  27Ti  Ax  $x(a-a) 

j  =  -mx  l-e  * 

Sin7T  AXSX  (2  mx*l)(a  -5) 

Sin7T  Ax  Sx(a-a) 


This  formula  holds  for  a  /  a  .  For  a  =  a  the  surn  is  evidently 

(2m  +1). 

x  1 

Now  if  we  take  A  6  =  - -%  the  numerator  is  zero  for 

xx  i  m  +1 

all  the  (integral)  values  of  (a  -  a  )  while  the  denominator  is  not.  From 
this 

27Ti  [AXSX  j(a-a)+  AySy  k  (b-b)  j 

j  k 


I  Z  e 


(4.  1) 


/■  (2mx*l)|2my*l)  if  a=a,b=b 


V.  0 


otherwise 
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This  means  we  can  solve  (4.1)  for  the  c  ,  in  terms  of  the 

a,D 

c  j.k  when 

- 1 -  - - £ -  l 

(2mxM)Ax  (  2  my  + 1)  Ay  J 


Let  us  display  the  two  formulas 


m_x  mx 


C',k  =  k?  C°>b  6 

o  —  —  b  =  - 1 1 1 y 


o-rr-  (  01  bk  \ 

<Lir'  \  2mx+l  ’  2my*l  / 


(4.2) 


'I 


/•  I  V'  w  ^  -27Ti(-~°1  ■■  +  -  bk  -) 

'’•b'  (!vO(!Yl)  ?£Ci,ke  * 


(4.3) 


°  12^.1)  (Z-r'Tn  Hci.kCos27r(T^7+  1^7?) 


These  relations  are  analogous  to  the  usual  reciprocal  relations  for 
Fourier  series,  e.g.,  (2.5)  and  (2.8).  More  important  to  the  prob¬ 
lem  of  filtering  is  the  following  convolution  formula.  Suppose 
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C  and  C  are  given  by  c  and  c  through  formula  (4. 2),  then 
their  convolution  is  defined  as 


rrv 


(  2mx  +  I )  ( 2my  +  l ) 


’x 

Z 

a=-mx 


17V 


I  I  ly 

1 

/3  =  -my 


Ci-a,k-/3  ca/3 


=  T; 


jk 


(4.4) 


where  C  is  regarded  here  as  a  periodic  function  having  values 
outside  the  range  (-mx.fflx)  To  be  specific,  Cmx  +  l,k' 

c-mx,k  •  cmx  +  2,k  =  C-mx  +l,k;etc.  If  we  transform  Tju.  by  a 
formula  such  as  (4.  3) 


b‘  (2mx  +  l)2(2my.  I  )*  Ci-ak-sCa^  e 


ot,(5»L-+ 

\  2  m_+  I  2m„+  I  / 


- 5 - 2  HC  e  277"'  ^  2"V  +  I  +  2"V'  )* 

(  2mx  +  I  )z  (  2my  +  I  )2  £  f  LCLf3  e  y 


z 

k 


Cj-a  k  -  /3 


e 


27Ti 


( q(i-ct) 

'  2mx+  I 


+ 


M* =J3±) 
2my+  I  / 


Since  in  the  second  sum  both  "C"  and  the  exponential  term  are 
periodic  with  period  (  2mx+l  )  by  (  2my+l  )  we  can  rearrange  the 
summands  toachievethe  sum  that  could  be  gotten  by  erasing  a  and  /3  . 
This  gives  us,  from  (4.  3), 


cab  cab 


(4.5) 
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So  th'j  Fourier  coefficients  of  a  convolution  are  the  products  of  the 
coefficients  of  the  two  functions  being  convolved. 

If  we  define  a  convolution  of  coefficients  as 


ca/3  =  yob 


(4.6) 


then  we  may  derive  similarly, 


I  I  w 


27Ti  ( 


a  J 


bk 


2mx  ♦! 


2  my+ 1 


+  1  )  5  Cik  Ci 


jk  ^jk 


jk 


(4.7) 


so  that  the  expansion  whose  coefficients  are  convolutions  is  the 
product  of  two  expansions. 

We  will  apply  (4.4)  and  (4.5)  to  the  situation  where  ca  b  = 

r  ,  of  (3.1)  so  that  c.  ,  are  spectral  estimates  of  the  type  (2.2) 
a ,  b  J  t 

whose  filter  function  is  (2.  6).  With  this  we  may  take  c  ,  so 


'a  b 


COS  27 r ( 


aj 


bk 


2mx*l 


2  my  ♦  I 


)=w(a,  b  ;  j  Sx  ,  kSy) 


(4.9) 
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* 


where  w  are  the  coefficients  for  the  spectral  estimates  (3.3),  and 
correspond  to  a  desirable  filter  function  W('),  (2.5).  Then  from 
(4.  4)  and  (4.5)  and  (3.3)  we  have, 


a,0 


E  E  w(a,b;a  8X  .^Sy)rab  =  E  E  cq5cos277-(' 
a  b  a  d 


aa 

2mx+l 


b  £ 


+?)' 


2my+l  /  a  b 


(4.  10) 


2mx  + 1 


2m, 


’  + 1  ?  ?  Ca  -  a  ,  /3  - 


(3-b  a  b 


From  this  we  have  a  fundamental  relation.  C  is  the  "raw  spectrum" 

of  the  covariances  r  .  and  the  coefficients  C  arc  used  to 

a,b 

"smooth"  it  by  taking  a  running  average.  By  this  we  may  achieve 
a  variety  of  W(‘)  functions.  Note  that  the  actual  W(')  function 
is  a  transform  (2.5)  of  the  product  of  the  C's  with  cosines  and  is 
therefore  a  convolution  of  the  type  in  (4.8).  All  the  filters  used  in 
practice  are  given  either  in  terms  of  the  running  average  weights 
C  or  in  terms  of  the  coefficients  c  which  modify  the  covariances 
before  making  a  Fourier  transformation.  Corresponding  to  the  two 
sides  of  the  convolution  theroem,  we  have  two  ways  of  computing 
the  spectrum.  We  may  first  compute  a  "raw  spectrum"  or  cosine 
transformation  of  r&  ^  and  then  smooth  it  by  a  running  average. 

This  is  convenient  when  only  a  few  of  the  running  average  coeffi¬ 
cients  C.  ,  are  non-zero.  Or  we  may  modify  the  r  .  by  mul- 

tiplying  them  by  the  factors  c  .  then  take  the  cosine  transformation. 

a ,  u 

This  is  done  when  the  c's  are  easy  to  compute. 
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An  application  of  the  convolution  relations  is  to  smoothing 
the  profile  which  was  discussed  for  the  one-dimensional  case  in 
Section  II,  Part  4.  We  should  note  that  the  work  of  the  section  deals 
primarily  with  the  data  and  not  with  the  mathematical  model  of 
randomness.  The  various  sequences  r^,  f^,  are  handled  as 
even,  periodic  sequences.  In  actual  computation  of  the  sums  we  do 
not  go  back  to  the  beginning  after  exhausting  the  table  so  that  the 
numerical  results  are  affected  by  "end  effects"  which  cause  them  to 
depart  somewhat  from  the  formulas  presented  here. 

Section  5.  The  two-dimensional  spectral  computer  program. 

The  program  operates  using  data  as  outlined  in  the  pre¬ 
ceding  sections.  The  general  operation  of  the  program  involves 
the  steps: 

(1)  The  area  data  is  smoothed  by  taking  differences  from 
one  fitted  quadratic  polynomial. 

(2)  The  covariances  of  the  smoothed  data  are  computed. 

(.3)  A  preliminary  spectrum  of  these  covariances  is 

computed. 

(4)  A  trigonometric  polynomial  is  fitted  to  the  reciprocals 
of  the  spectral  value's,  its  coefficients  are  used  to  make  a  running 
average  smoothing  of  the  covariances. 

(5)  A  new  spectrum  of  the  smoothed  covariances  is  com¬ 
puted. 

(6)  The  spectrum  is  smoothed  by  the  Hamming  method. 

(7)  The  spectrum  is  corrected  using  the  reciprocals  of 
the  polynomial  (4)  as  factors. 

The  program  was  checked  by  running  some  sea  surface 
data  from  the  report  by  Pierson  mentioned  in  Section  I.  The  re¬ 
sults  checked  very  well  except  for  the  latter  part  of  the  program. 
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Recently  this  part  has  been  made  to  operate  correctly. 

The  model  of  this  program  was  the  one  (Section  II,  Part  5) 
which  we  were  using  on  the  linear  data  at  the  time  the  specifications 
were  sent  to  the  programmer.  Our  later  experience  with  the  linear 
program  indicates  that  the  quadratic  polynomial  method  of  smoothing 
is  not  enough  for  our  purposes,  and  that  the  prewhitening  steps  may 
be  accomplished  better  by  using  an  improved  method  of  smoothing.  In 
order  to  complete  the  two-dimensional  analysis,  the  program  is  being 
modified. 

6.  Preliminary  Spectral  Estimate 

Table  5  contains  a  preliminary  spectral  estimate  from  our 
present  program  using  the  Aberdeen  data.  These  results  are  repro¬ 
duced  to  illustrate  some  of  the  comments  we  have  made  regarding 
two-dimensional  results.  Complete  studies  of  this  two-dimensional 
data  will  be  presented  in  a  later  report. 
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TABLES  AND  CHARTS 

Each  table  of  results  consists  of  five  parts.  The  first  part  is  the 
raw  profile  heights  given  in  feet  and  hundredths  above  the  surveyors' 
datum.  The  heights  are  presented  in  reading  order,  left  to  right  and 
down.  The  numbers  of  feet  from  the  beginning  for  certain  of  the  heights 
in  the  first  column  are  given  on  the  left  margin.  On  the  top  margin  are 
given  the  numbers  of  feet  to  add  to  the  position  of  an  item  in  the  left  column 
so  as  to  obtain  the  positions  of  the  heights  in  other  columns  of  the  same 
row.  The  second  parts  of  the  tables,  labled  A,  are  indexed  similarly. 

These  parts  contain  smoothed  heights. 

The  third  tables,  labled  B,  give  autocorrelations  of  the  smoothed 
2 

data  in  ft.  for  the  various  lag  numbers,  k.  The  fourth  tables,  labled  C, 

2 

give  the  spectral  density  estimates  of  the  smoothed  data  in  ft.  for  vari¬ 
ous  frequencies,  X  ,  in  cycles/ft.  The  last  tables,  D,  give  the  corrected 
spectral  estimates,  and  are  labled  similarly. 

After  making  the  calculations,  several  small  errors  were  found  in 
the  raw  data  of  Aberdeen,  Knox,  and  Yuma.  The  number  of  the  errors 
was  small,  at  most  fifteen  in  the  Yuma  data,  and  the  sizes  were  on  the 
order  of  a  few  tenths  of  a  foot.  The  effects  on  the  calculations  were  deemed 
to  be  minor  so  they  were  not  done  over.  The  raw  data  portion  of  each 
table  is  presented  in  a  correct  version,  but  the  other  parts  are  as  calcu¬ 
lated. 
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TABLE  1C 

POWER  SPECTRAL  DENSITY  ABERDEEN 
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TABLE  2B 


AUTOCORRQLATION  KNOX 


k 

1-1900  ft. 

1-950  ft. 

950-1900  ft. 

_ y 

0.0298 

0.0313 

’ 

0.0282 

i 

0.0168 

0.0176 

0.0159 

2 

.  .  0*0062 

0.0063 

0.0062 

3 

- 

0.0036 

_ 

0.0027 

- 

0.0046 

A 

- 

0.0106 

__ 

0.0108 

— 

0.0103 

0 

- 

0.012A 

— 

0.01A3 

- 

0.010b 

6 

- 

0.0126 

_ 

0.01A2 

— 

0.01  1 1 

f 

- 

0.0103 

_ 

0.0122 

- 

0.008b 

H 

- 

0.0061 

_ 

0.0074 

- 

0.0048 

9 

- 

0.0009 

— 

0.0022 

0.000b 

_ 10  _ 

0.0028 

0.0040 

0.00  18 

1  1 

0.0063 

0.0084 

0.004b 

12 

0.0062 

0.0081 

0 . 0  0  A  A 

1  3 

0.006b 

0.0079 

0.0053 

.  1A _ 

0.0061 

0.0082 

O.OOAO 

lb 

0.0032 

0.00b9 

0.0005 

..  . 

_  0.0000 

0 . 00 1 A 

- 

0  .  C  0  l  A 

17 

- 

0.0033 

_ 

0.0027 

- 

O.OOAO 

ia 

— 

0 . OuA  9 

0.0070 

0.0029 

10 

- 

0.005A 

— 

0.0091 

- 

0.001 7 

. . 20 _ 

- 

0  •  0  0  A  7 

_ 

0.0087 

— 

o.ooo  r 

21 

- 

o.ooHb 

0.0076 

0.0007 

„  22  . 

- _ 

. . .0.0011 

— 

_ 0.0037 _ 

_ _ 

0.0018 

23 

0.001C 

0.0002 

C . 0020 

2b 

0.0032 

0.0050 

0.0016 

2b 

O.OObl 

0.0092 

0.0012 

_ 20..... 

_  0.0047 

0.0103 

0.0009 

2  7 

0.00 3 A 

0.0078 

0.0012 

_ 26  ... 

0.0012 

0.0037 

0.001 7 

20 

— 

(i.OOO  1 

0.0013 

- 

0.0021 

30 

— 

<  .0023 

_ 

0.0032 

- 

0.0020 

31 

— 

O.OOA3 

_ 

0.0071 

- 

0.0021 

32 

- 

_  O.OOAb. 

— 

_ 0.0092 _ 

_ 

0. 0000 

33 

- 

0.0033 

_ 

0.0085 

0.0021 

3  A 

- 

. .  0.0011 

_ 

. 0,0.055 _ 

0.0039 

3  b 

0.000  3 

_ 

0.0018 

0.0030 

36 

0.0021 

0.0029 

0.J019 

37 

0.0037 

0.0067 

0  .  U  j  1  A 

3a 

0.00  3  A 

0.0067  .  .  . 

_ 

0.0005 

<0 

('.('030 

0.0073 

— 

.0.0  009 

AO 

C  .002  l 

0.0071 

- 

0.0029 

A  1 

O.UUOA 

_  0.0039 _ 

0 . 0  )  3  5 

A2 

- 

0.0018 

- 

o.oo  u 

0.0033 

A  ^ 

_ 

0,0026.  _ 

— 

0.0049 

_  “ 

0.0015 

A  A 

- 

0.0029 

- 

.  0. 0076 

0.0006 

Ab 

- 

0,00  2  7 

- 

0.U081 

0 .00  la 

A6 

- 

0.0019 

- 

0.0071 

. 

0.0029 

A  7 

- 

it.  00  13 

- 

0.0050 

0 . 002  5 

AH 

- 

().C")0  7 

- 

O.uulf 

0 .0007 

AO 

o.OnOb 

0.0017 

0.0000 

bO 

0 . 002 1 

0.00  72 

- 

0.9019 
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TABLE  2C 


POWER  SPECTRAL  DENSITY  KNOX 

SMOOTHED  DATA 


X 

-1-1900  ft. 

1-750  ft. 

750-1900  ft 

-0 .00000 

0.C0042 

0.00671 

G.CU370 

,  0.005Q0 

0.00720 

0.01029 

0.00320 

0.01000 

0.01290 

0.C0370 

0.0191 2 

.  0.01500 

0 . 04368 

0.C4129 

0.04414 

0.02000 

0.07294 

0.05083 

C.C8900 

_ 0 .02  5_00_ 

0.  1644  3 

0. 10686 

0.20543 

0.03000 

0. 30H29 

0.25217 

0.36722 

..0,0  3500. 

0. 70377 

0.91519 

0.4  /9 2 3 

0.04000 

0.95972 

1 .23784 

0.  7 0 3 0 f. 

..0.04  500 

..0.62  756 

0.56708 

0.70781 

0.0*3000 

0.37234 

0.31496 

0.42663 

_ 0. 0  5  5.00  _ 

0. 32066 

0.34175 

0.25879 

0.06000 

U. 25 758 

0. 29263 

0. 22290 

0 .06  500 

0.20449 

0.19727 

0.2093  1 

0.U7000 

0.  1  7  2  24 

0.  13977 

0.20400 

__0 .0  7*>00  _ 

0. 14001 

0.098.33 

0.18329 

0.0H000 

0.11015 

0.C9414 

0. 12397 

_ 0 .08  5.0  0_ 

0 . C 9245 

0.C9207. 

0.09176 

0. 09000 

0.08590 

0.08395 

0.08302 

„ 0.09500 

_ .0 ...09092 

0. 11834 

-0...C  7>3  3.4 

0.10000 

0.  10691 

6.  12  719 

C. 08727 

...0,10500 

0.08625 

0.06292. 

_0.  109  2.3_ 

0.1  1000 

0.06753 

0.0  3  763 

0 . 0  7  6 1>  4 

_0. 1  15Q0__ 

0.04  4  n 

0.C3  153 

JLDJU  2JL 

0.1  2000 

0.04182 

0.03970 

o .  o  4 1  n 

..0-12500 

0.02  7o7 

.0.0  31  19 

0.02227 

0.1  1000 

u.C 3201 

0.04644 

0.01699 

..0.1  3500 

9.04  6  3  5 

O.C  7  903 

0.0 1 7  ;2 

0.14000 

U . 0  4  166 

0.0  7  103 

O.C 1261 

__0 . 14500 

O.C  1002 

.0  .  C  4  9  ?  7_ 

0.0 1034 

o.isooo 

0 . 0  3  2  63 

0. C  5482 

0.01078 

0.15500 

. 0.02  7 32 
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TABLE  3B 

AUTOCORRELATION  YUMA 


k  1 -1500  ft. _ 

o  _ o.'.n*2 

-r - -  o .  o  1 1  o 

?  0.0031 

,  _  0.003d 

,  _  0.0096 

\  _  0.0040 

_ _ _ o  .0027 

— j  I  u .  0  0  0  I 

'  0.0004 

-  0.00 l 3 

0  0.UU14 

,  -  0.0002 

\\  _  0.0001 

— *  , - 1  0 . 0  0  L  o 

_  o.oo io 

U - 0.0014 

0.0U04 

},-■  .  _  0.0006 

0 .000(1 

— - 1  '  0.00  14 

_  0.00  10 

..  '"1"'  0.0005 

0.0002 

",  0.0015 

'u  0.0014 

—ti - :  o.oooi 

...  _  0.0016 

y.  _  0.0022 

_  0.0023 

, ;  -  0.0(  14 

v,  -  o.yooi 

- ,  0.0005 

‘  0.00 H 

0.0014 
0.0006 

\l  -  0.0004 

_ _ o.oo  iO 

■ - -  __  0.0003 

_  0.0007 

. .  . 

>  n  0.UOC4 

0.0005 
0.0000 

y,  -  o.oooi 

’  0.0004 

J:  o.ooo? 

A  O.OOOB 

4  6  -  °-U0^ 

;  -  O.LOO't 

..  _  0« 

;  -  o.ooo  i 

-U -  0.0004 


1-750  ft. 

0.0109 
0.0104 
(1 . 1)06 1 
O.OuOO 
0.0026 
0  .  u  0  3  6 
0 . 0042 
0 . 1)02  3 
0.0016 
0.0001 
O.O0C9 
0.00  1  7 

(.0017  _ 

0.0007 
0.00  1  l 
0. J030 
0 .0021 
0.0013 
0.0012 
f  / .  0  C  0  4 
0.0001 
C .0004 
0.0  UO  7 
o.OOOd 
0.00  15 
0.002  3 
Ci.i.024 
o . 0020 
O.001U 
r.00  1 1 
0 ,0002 
O.0003 
0.000  1 
0.0002 
0.0001 
0.0001 
0.0003 
0.0009 

_ 0.0  OJ? _ 

0.0012  — 
0.0014 
C . 0005 
O.oOOO 
0.0CQ2- 
r . oo  1 2 
0.0012 
C  .00  20 

0.0913 

r.ooi  i 

0.0004 

0.0004 


750-1500  ft. 

_ 0,0  24  5 _ _ 

0.0117 

0.0005 

0.0004 

0.Q0H6 

0.0044 

_ Q..OOJ.R - 

0.0009 
0.0024 
0.0025 
0.0019 
0.0013  1 

0 .  UO  2  3 _ J 

o.oo? a 
0.0009 
0.0001 
0,0011 
0.0004 

_ _Q,.QQ04 - 

0.0026 

ofoo3a.  - 

0.0014 

.0.0013 

0.0039 

_ _ 0,210.45 - 

0.0021 
0.0009 
0.0025 
0 . 0030 
0.0019 

u .  0.Q.05 - 

0.0007 
0..0029  .  - 
0.0032 
0,0013  — 

0 .000(1 

_ 0.1IQ2i - 

0.0021 

0^0Q26._ 

0.0QQ5- - 4 

0.0003  ( 

0.U0U5-  — 
0. 1)002 

_ U.QULO - 

0.0014 

0.UQ2U . 

u.ocon 

0 . 0004  .  . 

0.0005 

_ 0.1. OOP  -- 

0.0007 
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TABLE  3C 


POWER  SPECTRAL  DENSITY 

YUMA 

X 

1-1500  ft. 

1-750  ft. 

750-150  ft 

-0. 

C.C3108 

0.C3456 

0.62306 

..  0.00500 

O.C 7247 

0.  109  78 

0..  0.3.3  /? 

0.0100C 

0.13516 

0.20851 

0.06286 

0.0150C 

0.  14861 

0.20156 

0. C9640 

0.02000 

C.  13675 

0.  14153 

0. 13135 

0.02500 

0.  14893 

0.19571 

0.1013? 

0.03000 

0.  15961 

0.22116 

0. 09656 

Si  0.0  3500 

0.17715 

0.22174 

0.12865 

'  0.04000 

0.24567 

0.2/879 

0.21640 

0.04500 

0.32252 

0. 33020 

0. 32251 

0.05000 

0.27873 

0.7174? 

0.241/5 

0 .05500 

_0_.  2  1.49.5. 

(i._189.69_ 

0. 27.628 

0 .06000 

0.26751 

0.  10  34  3 

0.  4 7 3 8 ° 

...0 . 0650C 

...0.26  7  37. 

..0...C3 /I  l. 

0.45017 

0.0/00!; 

0.  16216 

0  .  C  .3  2  6  1 

0. 7413? 

__0 .07500 

.0.  119U8 

_0.  C8064. 

0.  15.6  7  4 

o.oaoco 

0.13259 

0.090^2 

C. 1 7507 

_ 0 .08500  _ 

0. 1 353? 

0.0 7840 

U.  19 32 9 

0.09000 

0.  12021 

0.05434 

0 .  18  5  7  8 

__  0. 09  500 

.0.0  )  1  >2. 

0_.  0  36  37 

.0.9  2.5  >6 

0.10000 

0.0/ 369 

0.0288/ 

C. 1  lb  15 

...0.10500 

..  0.0  / 26  3. 

O.C?  20  6. 

o .  1?  3?  3 

0.1 1000 

O.Ou45l 

0.C2919 

( .0994 7 

_0. 11500  _ 

_.0 .03?  5.9. 

J0..04  39JL 

.'  .06104 

0.12000 

0.04383 

0.04238 

O.C'-4  3  7 

._  0.12500  . 

1).  0  3  369 

0.0 3544 

C. 031 17 

0.1 3000 

0.0/4 >7 

O.C  3C30 

01783 

0.13500 

0.0._>  34  3 

0.01989 

(.('2712 

0.14000 

0.024 >6 

0.01864 

O.C? 1 49 

„  0. 14500 

O.C  2  6  8  7. 

.0 . 0  3?4«_ 

0.02077 

0. 15000 

0.0  >062 

0.039/9 

O.C,,?  1  39 

..  0.15500 

0.0279 > 

0.0260  7. 

0.01.969, 

0. 16000 

C. 0226 3 

0.02317 

0.02124 

0.16500 

O.C  2  >'-75 

O.C? 5  36 

O.i  3  2  3? 

o. i /ooo 

0.0 2 9u5 

0.C2 109 

0.C252! 

0.1/50i; 

0.01 7  2  < 

0. 02008 

0.01447 

0.  uoco 

0.0?  1  12 

0.02507 

0.0164! 

0.1*1509 

0.0  3 1 32 

"'•03  349. 

9.02931 

.  0.1900" 

0 .03/./ 

.03681 

3  861 

0. 1  >60'. 

O.C?/  JO 

9.02959 

0.02741 

0.20001- 

C.Cj'lfi'J 

•.'.02  7  59 

6.0  3 '3  4! 

0.20500 

0.9  3  >66 

0.02 >45 

3.03741 

_  0.2100  v<  . 

0  .  2  1 9. 

0.0 ->606. 

0.03(1*  7 

0.2 1500 

().:.••  /  /4 

0.0537? 

0.0421? 

0.22000 

0  .  O'.  09; 9 

0 . 04846 

0.03335 

0. 22500 

C. 42683 

0.03541 

o .  0  1  7  j  l 

0.23000 

0.9?  14  7 

0.02730 

O.C  18  6' 9 

0.23500 

0.O2W5 

0.02803 

0.02322 

0.24000 

".(,  ;0u4 

0.02.5/1. 

".03422 

0.2450" 

O.C  >  » J  H 

".02217 

0.04  8  3C 

0.25000 

0.0  it,:,  3 

O.C  1*)  36 

C.o4  3;0' 
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TABLE  3D 
POWER  SPECTRAL  DENSITY 


A 

0-750  ft 

0.  000 

— 

0.  005 

418. 09392 

0.  010 

50. 31736 

0.  015 

11. 36420 

0.  020 

3. 49368 

0.  025 

1. 67188 

0.  030 

. 94115 

0.  035 

. 63102 

0.  040 

. 57140 

0.  045 

. 53538 

0.  050 

. 35276 

0.  055 

. 23986 

0.  060 

. 24876 

0.  065 

. 20486 

0.  070 

. 11446 

0.  075 

. 08025 

0.  080 

. 08813 

0.  085 

. 09131 

0.  090 

. 08448 

0.  095 

. 06128 

0.  100 

. 05968 

0.  105 

. 06446 

0.  110 

. 06312 

0.  115 

. 05666 

0.  120 

.  05158 

0.  125 

. 04263 

0.  130 

. 03272 

0.  135 

. 03229 

0.  140 

. 03405 

0.  145 

. 03568 

0.  150 

. 03875 

0.  155 

. 02730 

0.  160 

. 02499 

0.  165 

. 02946 

0.  170 

. 02199 

0.  175 

. 01548 

0.  180 

. 01814 

0.  18  j 

. 02585 

0.  190 

. 03009 

0.  195 

. 02327 

0.  200 

. 02494 

0.  205 

. 02803 

0.  210 

. 03663 

0.  215 

. 04388 

0.  220 

. 03975 

0.  225 

. 02779 

0.  230 

. 02588 

0.  235 

. 03001 

0.  240 

. 03665 

0  245 

. 04396 

0.  250 

. 03901 

0-375  ft 

YUMA 

376-750  ft 

— 

.... 

633. 34277 

194. 53742 

77. 62410 

23. 40152 

15.  41482 

7. 37170 

3. 61580 

3. 35572 

2.  19704 

1. 13741 

1. 30409 

. 56937 

. 78989 

. 45790 

. 64843 

. 50332 

. 54813 

. 53536 

.  40172 

. 30621 

.  19356 

. 28205 

. 08949 

. 41002 

. 06674 

. 34492 

. 05831 

. 17039 

. 05436 

. 10564 

. 06016 

. 11636 

.  05290 

. 13043 

.  03854 

. 13056 

.  02717 

. 09422 

. 02338 

. 09570 

.  01958 

. 10965 

.  02856 

. 09732 

. 04731 

. 06578 

. 05046 

. 05221 

. 04485 

. 03939 

.  04036 

. 02375 

. 02717 

. 03705 

.  02546 

. 04296 

.  04380 

. 02758 

. 05024 

. 02707 

. 03095 

. 02338 

.  02558 

. 02345 

.02598 

. 03311 

.  02012 

. 02404 

.  01798 

. 01292 

. 02133 

. 01402 

. 02742 

. 02403 

.  02956 

. 03100 

. 02367 

. 02232 

. 02234 

. 02752 

.  02494 

. 03116 

.  03999 

. 03375 

.  04917 

. 03855 

.  04711 

. 03242 

.  03668 

.  01855 

.  03010 

. 02083 

.  03267 

. 02706 

.  03132 

. 04169 

.  02778 

. 06053 

.  02323 

. 05552 
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TABLE  4 

BATTLEFIELD  DAY  -  6201  FT. 

0  1  2  3  4  5 

ONE  FOOT  SPACING 

6  7  8 

9 

10 

11 

00003 

00009 

00016 

00022 

00027 

00035 

00042 

00048 

00053 

00059 

00064 

00069 

00074 

00079 

0008  1 

00084 

00098 

00  109 

00  116 

00  1  22 

001  30 

00139 

00146 

00152 

00165 

00174 

00180 

00189 

00195 

00201 

00206 

002  14 

00224 

00223 

00239 

00245 

00252 

00253 

00255 

00264 

002  7  1 

00277 

00286 

00294 

0030? 

00309 

00319 

00326 

00331 

00335 

00344 

00350 

00356 

0037  1 

00394 

004  13 

00429 

00433 

00439 

00445 

00459 

00478 

00491 

00498 

0051  3 

00527 

00534 

00540 

00549 

00557 

0056e 

00578 

00589 

00597 

00605 

006  1  8 

00635 

00648 

00664 

00682 

00700 

00709 

007  19 

00730 

‘00736 

00737 

00728 

007  !  8 

007  1  1 

0070  1 

00692 

00687 

00683 

00678 

00674 

00671 

00666 

00660 

00657 

00640 

00642 

00638 

00637 

00636 

00637 

00630 

00626 

00623 

00617 

00610 

00606 

00602 

00596 

00591 

00586 

00580 

00576 

00572 

00564 

00561 

00546 

00543 

00539 

00528 

0051  7 

00510 

00509 

00512 

00518 

00517 

00518 

00519 

00520 

00520 

00518 

00512 

00512 

00510 

00507 

00506 

00505 

00540 

00513 

00527 

00533 

00545 

00560 

00573 

00588 

00608 

00634 

00653 

00666 

00669 

00660 

00651 

‘00640 

00627 

00613 

00596 

00572 

00557 

00551 

00548 

00540 

00535 

00535 

00533 

00526 

0051  7 

00502 

00484 

00468 

00451 

00432 

00409 

00390 

00376 

00367 

00347 

00339 

00  332 

00327 

00324 

00324 

00324 

00327 

00324 

00320 

00317 

00314 

00308 

00307 

00300 

00295 

00292 

00287 

00283 

00279 

00275 

00273 

00268 

00266 

00268 

00269 

00275 

00281 

00285 

00293 

00296 

00296 

00296 

00294 

00299 

00303 

00302 

298 

00294 

00290 

00288 

00288 

00288 

00286 

00286 

00290 

00291 

00295 

00297 

00297 

'00298 

00300 

00302 

00305 

00303 

00300 

00296 

00297 

00298 

00297 

00296 

00297 

00301 

00302 

00302 

00302 

00302 

00304 

00307 

00309 

00309 

0031  1 

00315 

00316 

00320 

00323 

00324 

00324 

00327 

00330 

00332 

00335 

00335 

00334 

00331 

00330 

00328 

00326 

00326 

00325 

00335 

00335 

00335 

00334 

00336 

00337 

00339 

0034  1 

00340 

00334 

00332 

00336 

00340 

00342 

00343 

00346 

00340 

00336 

00340 

0034  1 

00341 

00340 

00338 

00336 

00335 

00336 

00338 

00338 

00344 

00346 

00344 

00343 

‘00345 

00349 

00353 

00349 

0035G 

00351 

00352 

00355 

00360 

00364 

00363 

00375 

00388 

00398 

00404 

004  1  0 

004  1  5 

004  1  8 

00421 

00426 

00428 

00426 

00422 

00416 

004  1  4 

004  1  1 

0041  1 

00412 

00402 

00393 

00305 

00381 

00370 

00374 

00370 

00368 

336 

00365 

00365 

00367 

00368 

0037  1 

00374 

00376 

00381 

00385 

00391 

00398 

00403 

00406 

00406 

00403 

00401 

00400 

00  395 

00388 

00381 

00375 

00371 

00365 

00363 

00364 

00363 

00364 

00364 

00367 

00374 

00301 

00386 

00391 

00390 

00391 

00389 

*00392 

00391 

00388 

00381 

00377 

00375 

00378 

00377 

00381 

00386 

00393 

0040  1 

00419 

00432 

00443 

00453 

00467 

00480 

00482 

00479 

00476 

00475 

00470 

00460 

00451 

00442 

00434 

00432 

00432 

00450 

00445 

00445 

00446 

00440 

00435 

00433 

00431 

00430 

00428 

00431 

00432 

00433 

00437 

00448 

00457 

00469 

00477 

00490 

00505 

00505 

00510 

0051  1 

00513 

00516 

00523 

00527 

00528 

00530 

00537 

00540 

00541 

00551 

00567 

00572 

00582 

00585 

00587 

00588 

00591 

00597 

00600 

00603 

456 

"00609 

00616 

00622 

00628 

00633 

00636 

00636 

00646 

00659 

00668 

00678 

00688 

00694 

00700 

00706 

00712 

00723 

00732 

00727 

007  17 

00712 

00709 

00707 

00705 

00709 

00713 

00717 

00721 

00722 

00727 

00725 

00726 

00728 

00733 

00740 

00743 

00747 

00752 

00757 

00758 

00  759 

00761 

00759 

00  756 

00758 

00763 

00769 

00779 

00783 

00791 

00801 

0081  1 

00822 

00833 

00842 

00849 

00856 

00857 

00854 

00850 

00845 

00834 

00823 

00816 

00816 

008  1  3 

00812 

00812 

00814 

008  1  6 

00820 

00825 

”00825 

00826 

00826 

0082  1 

00816 

00810 

00804 

00804 

00800 

00820 

008  1  7 

00817 

00021 

00834 

00840 

0083e 

00835 

0084  1 

00844 

00851 

00860 

00868 

00876 

00883 

00895 

00909 

0091  7 

0092  1 

00921 

00920 

00913 

00898 

00877 

00858 

0084  1 

00832 

00831 

00832 

00830 

00836 

00846 

00858 

00869 

00885 

00898 

0090  1 

00902 

0089  r 

564 

00889 

0088  1 

00863 

00848 

00838 

00832 

00829 

00828 

00828 

00833 

0084  1 

00844 

00845 

00846 

00844 

00850 

00855 

00854 

00856 

00860 

00863 

00866 

00866 

00866 

00073 

00879 

00881 

00883 

00885 

00887 

00807 

00883 

00872 

00865 

00860 

00857 

00852 

00846 

00843 

00842 

00848 

00851 

00857 

00869 

00879 

00887 

00894 

00897 

00898 

00900 

00898 

00892 

00887 

00884 

008  78 

00076 

00878 

00080 

00883 

00885 

00888 

00892 

00896 

00899 

00904 

00906 

00906 

00898 

00895 

00896 

00896 

00898 

00904 

0090  7 

00909 

00912 

009  1  4 

009  1  7 

009  1  8 

0092  1 

00923 

00923 

00920 

009  1  4 

00908 

00908 

00902 

00900 

00895 

009  1  0 

0092  1 

00922 

0092  1 

00926 

00930 

0U935 

00938 

00944 

00948 

00952 

00952 

0  0954 

00956 

00959 

00960 

00964 

00970 

00976 

00982 

00992 

0  1  001 

0  1  008 

01017 

0  l  025 

0  1030 

0  1031 

0  1  028 

01024 

01021 

01017 
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084 


804 


912 


1032 


1140 


1260 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

01013 

01011 

01012 

0  1  020 

0  1  030 

0  1  040 

01045 

01058 

01069 

01079 

01096 

01105 

01107 

0  1103 

0  1  099 

0  1  097 

01091 

01084 

0  1  074 

0  1  063 

01058 

0  1055 

0  1  053 

01053 

0  1052 

0  1  052 

01055 

61  058 

0  1  082 

Cl  069 

0  1  078 

01  0B9 

OT  106 

01130 

0  1149 

01163 

01176 

01185 

01190 

0  1192 

0  1194 

0  1191 

0  1189 

01185 

01180 

0  1179 

01181 

01187 

011  94 

01198 

01201 

01203 

0  1203 

01204 

01203 

01202 

01203 

01205 

01207 

01209 

01214 

01220 

01234 

0  1246 

0  1254 

01262 

01269 

01277 

01284 

0  1294 

0  1301 

01  337 

01312 

01315 

01318 

0  1  320 

01321 

0  1321 

0  1  320 

01320 

0  13  19 

01316 

01321 

01  323 

01332 

01  342 

0  1  350 

0  1  354 

0  1361 

01369 

0  1  372 

0  1377 

01379 

0  1381 

01383 

01  390 

01392 

01398 

01405 

0  1  409 

01417 

01426 

01437 

01  446 

01453 

01458 

0  1  464 

01467 

01466 

01465 

01465 

0  1  464 

01463 

0  1466 

0  1  470 

01476 

01479 

01487 

01507 

01522 

01535 

01580 

0  1562 

0  1572 

0  1581 

01582 

01585 

01588 

01586 

01583 

0  1578 

01579 

01580 

01587 

01596 

0  1603 

01606 

0  1602 

0  1593 

01591 

01585 

01584 

0  1585 

01588 

01593 

01601 

01611 

0  1  628 

01647 

0  1663 

01  680 

0  1693 

01700 

0  1705 

0  1709 

01713 

01714 

01712 

0  1  706 

0  1701 

0  1694 

0  1687 

01682 

01683 

01684 

0  1685 

0  1685 

01686 

01688 

01690 

01695 

0  1  700 

0  1  705 

0  1715 

01723 

01730 

01737 

01740 

01742 

0  1746 

0  1751 

0  1  754 

0  1755 

0  1  756 

01757 

0  1  760 

01760 

01759 

01761 

01762 

0  1765 

0  1769 

01771 

01772 

01770 

0  1768 

01789 

0  1  792 

01798 

01801 

01804 

01804 

0  1  806 

01814 

01821 

01824 

01832 

01847 

0  1  868 

0  1  893 

01911 

01925 

01939 

01955 

01963 

01962 

01  956 

0  1  946 

0  1  939 

0  1931 

01922 

0  1911 

01909 

01925 

01931 

01936 

01944 

01947 

0  1952 

0  1958 

01964 

01971 

01998 

02006 

02027 

02036 

02034 

02032 

02030 

02026 

02018 

02017 

02018 

02013 

0201  4 

02016 

02017 

02017 

02016 

020  1  3 

02008 

02010 

020  1  3 

02019 

02028 

02035 

02040 

02044 

02052 

02063 

02073 

02083 

02095 

021  03 

02110 

02117 

02126 

02  134 

02141 

02  1  44 

02  1  46 

02  1  48 

02149 

02155 

02160 

02165 

02173 

02  186 

02  196 

02207 

02218 

02227 

02237 

02241 

02245 

02249 

02257 

02265 

02274 

02285 

02294 

02305 

023  1  0 

02317 

02326 

02340 

02349 

02357 

02366 

02376 

02386 

02397 

024  1  2 

02423 

02432 

02444 

02458 

02468 

02477 

02487 

02497 

02505 

02512 

02520 

02526 

02536 

02551 

02565 

02587 

02609 

02629 

02645 

02647 

02644 

02639 

0263  1 

02621 

02606 

02593 

02585 

02576 

02569 

02562 

02560 

02559 

02562 

02567 

02575 

02588 

02591 

02593 

02596 

02598 

02602 

02605 

02605 

02606 

02607 

02617 

02624 

02634 

02647 

02662 

02684 

02707 

02722 

02732 

02744 

02752 

02753 

02  750 

02750 

02748 

02743 

0274  1 

02739 

02738 

02743 

02753 

02767 

02781 

02797 

0281  1 

0282  1 

02828 

02828 

0282  1 

02812 

02804 

02791 

02771 

02754 

02743 

02737 

02734 

02733 

02730 

02726 

02725 

02728 

02728 

02735 

02744 

02750 

02760 

02761 

02766 

02775 

02788 

02801 

028  1  7 

02834 

02845 

02853 

02859 

02863 

02866 

02867 

02865 

02862 

02860 

02859 

02861 

02862 

02859 

02856 

02863 

02865 

02870 

02881 

02886 

02893 

02904 

0291  1 

0292  1 

029  30 

02940 

02950 

02955 

02960 

02966 

02972 

02977 

02989 

03000 

0301  1 

0301  7 

03026 

03037 

03043 

03050 

03055 

03058 

03060 

03062 

03066 

03069 

03071 

03075 

03078 

03076 

03075 

03080 

03087 

0309  1 

03097 

031  1  1 

03123 

03131 

0314  1 

03153 

03165 

03173 

03  1  85 

03  1  89 

03196 

03199 

03209 

032  1  4 

032  1  4 

0322  1 

03225 

03231 

03236 

03242 

03251 

03258 

03264 

03271 

03277 

03282 

03292 

03294 

03300 

03305 

03310 

03309 

033  1  0 

033  14 

03317 

03325 

03334 

03347 

03360 

03372 

03377 

03384 

03399 

03409 

03416 

03423 

03427 

03431 

03436 

03445 

03446 

03452 

03455 

03459 

03464 

03465 

03467 

03468 

03469 

03474 

03478 

034  79 

0349  1 

03483 

034  76 

034  7  1 

03467 

03461 

03452 

03448 

03444 

03443 

03447 

03447 

03451 

03463 

03475 

03485 

03493 

03498 

03500 

03500 

03496 

0  349  1 

03487 

03484 

03482 

03480 

03477 

03480 

03483 

03487 

03495 

03504 

035  16 

03529 

03536 

03546 

0355  1 

03555 

0  3554 

03549 

0  354  0 

03535 

03535 

03533 

03532 

03532 

03530 

03529 

03527 

03529 

0  353  4 

03538 

03540 

0  354  3 

03553 

03561 

03568 

03581 

03599 

0360  1 

03606 

03610 

03609 

0  3610 

036  1  1 

036  1  3 

036  1  5 

03618 

03626 

03634 

03633 

03639 

03647 

03652 

03656 

0  3664 

03670 

03674 

03675 

03675 

03671 

03667 

03664 

03661 

03662 

03661 

0  3659 

0  3664 

03664 

03665 

03670 

03673 

03679 

03687 

03693 

03698 

03702 

03709 

0371  1 

0  371  3 

03713 

03712 

03715 

03713 

03713 

03714 

03718 

03723 

03730 

03735 

03738 

03740 

03746 

03  746 

03748 

03749 

03749 

03753 

03755 

03762 

03765 

03769 

03777 

03/87 

03792 

03794 

03794 

03798 

03800 

03802 

03805 

03807 

038  1  0 

03817 

03826 

0  3830 

03838 

03851 

03872 

03884 

03887 

03893 

03899 

03903 

03908 

03913 

0391  1 

03916 

03921 

03928 

03934 

03936 

03944 

03953 

03963 

03973 

03984 

03995 

04004 

040  1  6 

04029 

04036 

04039 

04039 

04043 

04052 

04054 
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1368 


1488 


1608 


1716 


1836 


1944 


0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

10 

11 

04059 

04067 

04066 

04093 

04112 

04  1  25 

04  1  34 

04141 

04  1  45 

04  1  56 

04  1  67 

04  1  76 

04  1  74 

04  1  74 

04  1  77 

04  179 

04  193 

04200 

04204 

042  12 

04217 

04226 

04246 

04251 

04258 

04271 

04281 

04287 

04288 

04292 

0430  1 

04309 

04315 

04323 

04330 

04340 

04352 

04360 

04371 

04381 

04  389 

04394 

04400 

04405 

044  1  5 

04425 

04434 

044  38 

04449 

04459 

04469 

04477 

04486 

04495 

04512 

04  524 

04535 

0454  1 

04551 

04562 

04575  -Q45B6 

04594 

04599 

04609 

0461  3 

04615 

046  18 

04628 

04635 

04643 

04653 

04658 

04669 

04685 

04693 

04702 

04709 

04716 

04728 

04737 

04748 

04752 

04766 

04776 

04787 

04802 

04B  1  5 

04B27 

04B4  1 

04852 

04864 

04878 

04B86 

0  4895 

04903 

04910 

04915 

04920 

04922 ' 

04928 

04937 

0494  1 

04953 

04957 

04966 

04977 

049B7 

04993 

0500  1 

05010 

05018 

05023 

05030 

05037 

0504  1 

05046 

05052 

05064 

05071 

05083 

05102 

05110 

05114 

05115 

05113 

05110 

05  1  07 

05104 

05095 

05086 

05082 

05083 

05091 

0510  1 

05118 

05132 

05  1  44 

05151 

05154 

05  1  54 

05151 

05145 

05137 

051  32 

05127 

05125 

05128 

05132 

051  38 

0  5  1  4  B 

05  1  59 

0516B 

05173 

05172 

05175 

05178 

0518  1 

05176 

05168 

05  164 

05  1  59 

05154 

05  1  54 

05152 

051  46 

05129 

05111 

05103 

0510  1 

05098 

05094 

05095 

05099 

05097 

05093 

O50B4 

05075 

05069 

05066 

05060 

05056 

05054 

05055 

05060 

05062 

05063 

05065 

05070 

05075 

05082 

05093 

05095 

05091 

05089 

05088 

05088 

050B5 

0508  1 

05076 

05073 

05069 

05064 

05057 

05055 

05050 

0504  1 

05035 

0502  7 

05014 

05002 

04992 

049B3 

0497  1 

04962 

04951 

04942 

04938 

04938 

04937 

04944 

04952 

04967 

04963 

04959 

04955 

04953 

04950 

04953 

04952 

0495  1 

04954 

04950 

04946 

04945 

04942 

04940 

04935 

0493B 

04940 

04926 

04923 

04918 

04914 

049  1  9 

049  1  2 

04909 

04907 

04916 

049  1  B 

049  1  4 

04914 

04910 

04906 

0490  1 

04905 

04903 

0490B 

04902 

04897 

04B93 

04BBB 

04BB2 

04B76 

0487  1 

04866 

04866 

04875 

04BB3 

04B98 

04908 

04909 

0491  1 

04906 

0490  1 

04910 

04918 

049  1  4 

04910 

04906 

04904 

0490  1 

04B94 

04893 

04887 

04883 

04879 

0487B 

04876 

048  7  3 

04872 

04874 

04BB1 

0490  1 

04916 

049  14 

0491  4 

049  1  6 

04914 

0491  1 

'  04915 

04930 

04927 

04933 

04934 

04935 

04939 

04949 

049B6 

04  9B  7 

04986 

04980 

04997 

04997 

04993 

04989 

04983 

04981 

04976 

04970 

04966 

04962 

04964 

04966 

04965 

04962 

04960 

04959 

04  95B 

04957 

04977 

04971 

04969 

04969 

049B3 

049B5 

04985 

04987 

04991 

04994 

04994 

04997 

04999 

04999 

0499B 

04993 

04990 

0498B 

04987 

04985 

04990 

04994 

04992 

05000 

05009 

050  15 

05023 

05024 

05020 

05021 

05020 

05020 

050  17 

050  1  0 

05007 

05004 

05003 

05007 

050  1  1 

050  1  5 

0  5019 

05029 

'  05035 

05044 

0  5052 

05053 

05053 

050AB 

05047 

05044 

05044 

05061 

05057 

05054 

05051 

05049 

05048 

05058 

05063 

05062 

05061 

05060 

05058 

05055 

0  5051 

05047 

05043 

05034 

0503  1 

05046 

05045 

05036 

05037 

05039 

05031 

05026 

0502B 

05030 

05027 

05032 

05044 

05048 

05049 

05068 

0508B 

05097 

05103 

05110 

05116 

05119 

05122 

05  120 

05  1  34 

0813  1 

05125 

05123 

051  29 

05118 

05  116 

05119 

05121 

05124 

05122 

05  1  19 

051  15 

05109 

05  1  09 

05  1  05 

05098 

05  1  05 

05108 

051  05 

05107 

05113 

"05125 

05  1  34 

05137 

05  1  54 

051  76 

0  5  1  96 

05203 

05218 

05226 

05227 

05235 

05244 

05253 

05258 

05262 

05265 

05268 

0  52  70 

05276 

05284 

05291 

05296 

05303 

053  1  0 

05324 

05345 

05356 

05  364 

0537  1 

053B2 

05394 

05405 

054  1  3 

05428 

05452 

05470 

05475 

05482 

0549  1 

055D3 

055  1  5 

05525 

05532 

05538 

05550 

05583 

05598 

05605 

056  14 

05635 

05648 

05655 

05658 

05670 

0568  1 

05690 

0570  1 

057  16 

05730 

05740 

05752 

05762 

0577  1 

05782 

05800 

05815 

05026 

05B39 

05848 

05858 

05B65 

0  5  B  6  7 

05869 

05886 

05889 

05086 

05680 

05878 

05870 

05863 

05857 

05054 

05859 

05869 

05893 

05899 

05904 

05916 

0592  3 

05927 

05932 

059  35 

05937 

05938 

0  5954 

0096  1 

05967 

05978 

05983 

05990 

08997 

0600  1 

06007 

06015 

0603  1 

06037 

06054 

0606  1 

06066 

06070 

06073 

06095 

06107 

06110 

06114 

06116 

06121 

06127 

06137 

06  1  45 

06  1  50 

06158 

06  167 

0  6  17  3 

06  1  76 

06  1  82 

06  1  89 

06  193 

06191 

06190 

06192 

06195 

06  1  99 

062  1  0 

0621  1 

062  1  7 

06240 

06247 

06253 

06260 

06271 

06279 

06304 

06308 

063  1  3 

0631  7 

06322 

06333 

06340 

06349 

06364 

06364 

0637  1 

0  63  76 

06382 

06392 

06399 

06**  0  4 

06400 

064  1  4 

064  1  4 

064  1  6 

064  22 

064  2  7 

06435 

06442 

06456 

06482 

0650  1 

065  1  1 

06515 

06520 

068  35 

06648 

06650 

06552 

0655  3 

06553 

06552 

06550 

06548 

0655  3 

0  6556 

068  7  0 

0657? 

06577 

06588 

06603 

066  1  4 

0  6620 

06627 

06634 

06638 

06653 

0666  3 

06668 

066B2 

0060? 

0  6  606 

06  70  1 

06706 

06708 

06709 

06709 

067  10 

067  1  0 

067  1  0 

06  7  1  2 

06  715 

ot.  r  i » 

06  7  70 

06  724 

067  3  1 

0  6  74  8 

06757 

06  765 

0677  1 

06784 

06  790 

06792 

Of.  79  3 

06  7<*>o 

0  6  7  OS 

06779 

0  67  75 

06778 

06778 

06  7Q1 

06762 

0r>7B  l 

On  78  3 

Of)  /tU't 

06  786 

06  79 7 

Of. 808 

06806 

06804 

06802 

06  795 

06  792 

06  786 

Of.  tn  i 

0  6  t  ft 

06  770 

06  7  79 

On  7  7»' 

067  7  7 

0  6  7  7  4 

06772 

06769 
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2064 


2184 


2304 


2436 


2556 


0 

1 

2 

3 

4 

5 

06768 

06766 

06759 

06754 

06748 

06746 

06739 

06737 

06734 

06730 

06724 

067  19 

06694 

06690 

06682 

06676 

066  72 

06666 

06626 

06622 

06620 

066  1  9 

06630 

06630 

06683 

06699 

06712 

06723 

06736 

06750 

068  1  8 

0683  1 

06847 

06856 

06  866 

06880 

06978 

06989 

06996 

07005 

07015 

07025 

07205 

07225 

07237 

07265 

07286 

0  7293 

07381 

07393 

07407 

0  7415 

07431 

07450 

0  74  86 

07488 

0  7507 

075  1  6 

07523 

07529 

076  1  8 

07626 

07632 

07637 

0  7643 

07651 

07685 

07685 

07694 

07704 

0  7711 

07727 

07758 

07759 

07762 

07765 

0  7  783 

07801 

07826 

07831 

0  7835 

0784  0 

0  784  1 

07846 

07913 

079  1  4 

0  7917 

07935 

0  7946 

07950 

07967 

07970 

07975 

07977 

07980 

07982 

08002 

08009 

08025 

08031 

08032 

08029 

08005 

08002 

0  7998 

08000 

08004 

08008 

08081 

08086 

08094 

08  1  05 

08119 

08125 

08204 

086  1  6 

08233 

C  8250 

0826  1 

08276 

08449 

08475 

08503 

08536 

08565 

08590 

08730 

08745 

08770 

08794 

08812 

08829 

08993 

090  1  3 

09038 

09072 

09105 

09129 

09322 

09347 

09368 

09390 

094  I  9 

09451 

09593 

096  1  2 

09638 

0965  1 

09667 

09688 

09836 

09856 

0  988  1 

0991  0 

09936 

00952 

10099 

10109 

10130 

10  147 

10  16  7 

1  0  1  74 

1  0259 

1  0269 

1  0279 

1  0295 

10314 

1  0324 

1  0488 

1  0508 

1  0540 

1  0569 

1  0568 

!  0  367 

1  0538 

10525 

10512 

1  0497 

1  0405 

10  4  7  1 

10372 

1  0365 

1  0358 

1  0350 

1  0344 

10340 

10326 

10320 

10313 

1  0304 

1  029 2 

1  0280 

10180 

10169 

10  159 

10  148 

10  132 

10123 

1  0008 

09993 

09980 

09966 

09949 

09938 

09856 

09047 

09828 

098  1  9 

09829 

09838 

09760 

09748 

09745 

09736 

09723 

097  1  2 

09665 

0966  1 

09664 

09657 

09650 

09643 

09592 

09583 

09578 

09570 

0  9559 

0954  9 

09498 

09488 

09495 

09502 

09505 

094  90 

09403 

09389 

09378 

09366 

0  9352 

09342 

09296 

09288 

09295 

09296 

0928  7 

09278 

09205 

09  1  95 

09187 

09  1  79 

0917  1 

09  1  64 

091  35 

09135 

09  135 

09  137 

091  JO 

09140 

09165 

09  169 

0918  1 

09183 

09197 

09  1  98 

09227 

09225 

09228 

09227 

0922  7 

0922  9 

09248 

09255 

09255 

09257 

0926  1 

0926  7 

0930  1 

09303 

09308 

093  1  0 

09318 

09  320 

0934  1 

09342 

09348 

0934  1 

09345 
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TABLE  4B 


AUTOCORRELATION  BATTLEFIELD  DAY 
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— , 
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NOTE: 

After  making  the  calculations,  several  small  errors  were  found 
in  the  raw  data  of  Aberdeen,  Knox  and  Yuma.  Since  their  effects  on 
the  results  would  be  minor,  the  calculations  were  not  done  over.  The 
raw  data  table  in  each  case  is  corrected,  but  the  smoothed  data  tables 
are  as  computed. 
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Figure  1  -  PROFILE  -  ABERDEEN 
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Figure  6  -  SPECTRAL  ESTIMATES 
FT.  KNOX 
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Figure  1 1  -  CUMULATIVE  DIAGRAM 
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